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I.  OVERVIEW 

This  research  program  has  focused  on  a  novel  approach  to  the  control  of  the  aerodynamic 
performance  of  lifting  surfaces  by  fluidic  modification  of  their  apparent  aerodynamic  shape,  or 
virtual  aerosurface  shaping.  This  flow  control  approach  emphasizes  fluidic  modificatiOT  of  the 
“apparent”  aerodynamic  shape  of  the  surface  with  the  objective  of  altering  or  prescribing  the 
streamwise  pressure  gradient.  Control  is  typically  effected  by  the  interactions  of  arrays  o 
synthetic  jet  actuators  with  the  cross  flow  that  displace  the  local  streamlines  near  the  surface  and 
thereby  induces  an  ‘apparent’  modification  of  the  flow  boundary.  The  operating  frequency  of  the 
control  jets  is  high  enough  so  that  the  actuation  period  is  at  least  an  order  of  magnitude  lower 
than  the  relevant  characteristic  time  scale  of  the  flow.  Therefore,  the  interaction  domains 
between  the  control  jets  and  the  cross  flow  are  quasi-steady  and  hence  the  induced  aerodynamic 
forces  are  virtually  time- invariant.  Earlier  work  at  Georgia  Tech  demonstrated  the  utility  of  this 
approach  to  flow  control  for  the  suppression  of  flow  separation  at  post  stall  angles  of  attack,  m 
contrast  to  other  approaches  for  control  of  flow  separation  that  rely  on  the  receptivity  of  the 
separating  shear  layer  to  time-periodic  actuation  and  on  Coanda-like  unsteady  reattachment  (e.g., 
Hsiao  et  al.  1990),  the  aerodynamic  forces  achieved  by  the  present  approach  are  virtually  tij^c- 
invariant.  More  recently,  the  utility  of  this  approach  was  demonstrated  in  the  absence  of  stall  at 
low  angles  of  attack.  For  example,  significant  reduction  in  form  drag  with  minimal  lift  penalty 
can  be  achieved  by  controlled  formation  of  trapped  vortices  through  the  interaction  of  a  synthetic 
jet  and  a  miniature  [0(0.0 Ic)]  passive  obstruction  near  the  leading  and  trailing  edges  of  lifting 
surfaces. 

The  present  report  is  comprised  of  four  sections.  Section  I  describes  the  interaction  between  a  jet 
actuator  and  the  cross  flow  over  the  surface  of  a  2-D  cylinder  with  the  objective  of  demonstrating 
the  alteration  of  the  streamwise  pressure  gradient.  This  section  also  compares  the  effects  of 
fluidic  interaction  with  that  of  a  surface-mounted  miniature  passive  obstruction.  Even  though  the 
details  of  the  interaction  domains  with  the  cross  flow  are  not  the  sarne,  a  properly-scaled 
obstruction  can  induce  the  same  global  aerodynamic  forces  as  an  active  fluidic  element  within  a 
limited  parameter  range  and  therefore  these  elements  can  be  used  to  complement  each  other.  The 
second  and  third  parts  are  concerned  with  the  modification  ofd  the  aerodynamic  perfonnance  of 
lifting  surfaces  at  small  angles  of  attack.  These  effects  are  explored  near  the  leading  edge 
(Section  II)  and  the  trailing  edge  (Section  HI).  While  near  the  leading  edge  the  effects  is  mainly 
the  formation  of  a  trapped  vortex  of  controllable  circulation,  near  the  trailing  edge,  the  actuation 
manipulates  the  Kutta  condition  by  a  miniature  Gumey-like  flap  integrated  with  a  jet  actuator. 
Finally,  Section  IV  focuses  on  an  investigation  of  the  fundamental  flow  mechanisms  of  the 
interaction  of  surface-mounted  synthetic  jets  with  a  turbulent  boundary  layer. 
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1.2.  THE  INTERACTION  OF  SYNTHETIC  JET  ACTUATORS  WITH  A  CROSS  FLOW 

Although  the  interest  in  control  of  separation  for  aerodynamic  applications  has  been  primarily 
focused  on  2-D  and  3-D  airfoils,  some  control  strategies  have  been  investigated  in  the  nominally 
two-dimensional  flow  around  circular  cylinders  (e.g.,  Williams  et  al.,  1991  and  Pal  &  Sinha, 
1997).  This  simple  geometry  is  particularly  attractive  since  it  minimizes  higher-order  effects  of  a 
specific  global  (e.g.,  airfoil)  geometry  and  has  the  distinct  advantage  that  the  baseline  flow  is 
extensively  documented  including  the  evolution  with  Reynolds  number  of  flow  separation  (e.g., 
Roshko  and  Fiszdon,  1969).  Amitay  et  al.  (1997  and  1998)  investigated  the  effect  of  a  surface 
mounted  synthetic  jet  at  Reo  up  to  130,000  over  a  range  of  jet  angles  and  showed  that  the  jets  can 
effect  an  increase  in  lift  of  up  to  Cz,  =  1  and  a  corresponding  reduction  in  pressure  drag  of  up  to 
30%.  The  structure  of  the  interaction  domain  between  a  jet  and  the  nominally  two-dimensional 
cross  flow  around  the  cylinder  were  investigated  in  detail  under  earlier  support  from  AFOSR. 
High-resolution  (0.133  mm)  digital  particle  image  velocimetry  (PfV)  was  used  to  investigate  the 
nominally  two-dimensional  (time-  and  phase-averaged)  flow  field  in  the  absence  and  presence  of 
actuation  within  a  domain  that  includes  the  synthetic  jet  actuator  and  the  separating  shear  layer 
on  the  top  surface  of  the  cylinder.  These  measurements  were  conducted  in  a  two-dimensional 
(0.91  X  0.051  m  cross  section),  closed  return  wind  tunnel  using  a  63  mm  cylinder  model 
equipped  with  a  2-D  jet  actuator  having  an  orifice  width  of  0.5  mm.  The  cylinder  can  be  rotated 
about  its  centerline  so  that  the  angle  between  the  jets  and  the  direction  of  the  free  stream  can  be 
continuously  varied. 

The  interaction  domain  between  the  jet  and  the  cross 
flow  at  high  actuation  frequency  is  inferred  from 
distributions  of  the  spanwise  vorticity  (Figures  I.la-c) 
that  are  measured  along  the  surface  of  the  cylinder  at 
Re£)  =  21,500  (to  provide  a  thicker  boundary  layer). 

The  time-averaged  vorticity  map  in  Figure  I. la  shows 
that  the  baseline  flow  separates  at  ^  =  90°.  Phase- 
averaged  vorticity  concentrations  in  the  presence  of 
actuation  (Figure  1. 1  b)  show  that  the  (counter  rotating)  g 
vortex  pairs  that  engender  the  jet  interact  with  the  wall 
boundary  layer  and  form  a  train  of  clockwise  (CW)  o.s 
vortices  that  are  advected  downstream,  become 
weaker  and  ultimately  disappear  within  three  to  four  o.4 
wavelengths  of  the  actuation  {X  =  O.SUJfc^.  These 
data  also  show  that  the  (upstream)  counterclockwise  Figure  I.l.  Normalized  vorticity  StDaci  =  4.0 

(CCW)  jet  vortex  is  accelerated  above  and  around  the  baseline  (a),  and  actuated:  phase-locked^,  and 

CW  vortex  and  rapidly  weakens  (within  one  time-averaged (c).(r=  63  ,c,  =  5.ixi(r ). 
wavelength).  The  time-averaged  vorticity  in  the  presence  of  actuation  (Figure  I.lc)  shows  a  finite 
interaction  domain  that  protrudes  into  the  cross  flow  above  the  edge  of  the  local  boundary  layer 
in  the  absence  of  actuation  and  ends  approximately  2-3 X  downstream  from  the  jet  orifice  and  is 
followed  by  a  thinner  boundary  layer  compared  to  the  baseline  fiow. 


x/D 
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The  interaction  domain  in  Figure  Lie  displaces  the  outer  flow  and  thereby  results  in  local 
variations  in  the  time-averaged  pressure  field,  which  is  computed  from  the  Reynolds-averaged 
Navier-Stokes  equations  using  the  PIV  data.  Figure  I.2a  shows  two  streamlines  having  the  same 
value  of  the  stream  function  in  the  absence  and  presence  of  actuation  that  are  superposed  on  the 
vorticity  distribution  of  the  actuated  flow  (note  the  displacement  of  the  streamline  in  the  presence 
of  the  actuation).  Here,  b  is  the  jet’s  orifice  width.  The  corresponding  pressure  distributions 
along  these  streamlines  are  shown  in  Figure  I.2b.  Compared  to  the  baseline  flow,  the  pressure 
along  the  streamline  in  the  presence  of  actuation  is  lower  over  the  entire  measurement  domain, 
indicating  an  increase  in  the  velocity  along  the  streamline  due  to  global  changes  in  the  flow  (e.g., 
the  movement  of  the  separation  point).  More  significantly,  within  the  azimuthal  sector  that 
includes  the  Interaction  domain  there  is  a  local  minimum  in  the  pressure  distribution  that 
indicates  a  significant  favorable  pressure  gradient  in  the  cross  flow  above  the  interaction  domain. 
The  favorable  gradient  is  followed  by  a  small  region  of  weaker  adverse  pressure  gradient,  and 
finally,  there  is  a  slight  favorable  pressure  gradient  through  the  rest  of  the  domain.  As  can  be 
inferred  from  Figure  I.2c,  the  alteration  of  the  streamwise  pressure  gradient  results  in  a  thinner, 
more  stable,  boundary  layer  downstream  of  the  interaction  domain  and  a  substantial  delay  in 
separation.  The  evolution  of  the  (time-averaged)  boundary  layer  displacement  thickness  S'  for 
the  baseline  and  actuated  flows  are  shown  in  Figure  I.2c.  Within  the  interaction  domain, 
(60°  <6<  80°),  6*  in  the  presence  of  actuation  increases,  reaches  a  local  maximum  immediately 
downstream  of  the  jet,  and  then  diminishes  substantially  suggesting  the  presence  of  a  favorable 
streamwise  pressure  gradient.  For  e>  80°,  d'  begins  to  increase  again  although  at  a  lower 
streamwise  rate  than  in  the  absence  of  actuation,  until  separation  occurs  at  110°. 


Figure  1.2.  (a)  Vorticity  in  interaction  domain  superimposed  with  streamlines  having  the  same  stream  function  value: 

baseline  (-—)  and  actuated  ( - );  The  pressure  coefficient  along  the  streamlines  (b)  and  the  boundary  layer 

displacement  thickness  (c):  baseline  (o),  actuated  (•). 


1.3  COMPARISON  WITH  A  PASSIVE  OBSTRUCTION 

Similar  to  the  effects  of  fluidic  actuation,  the  flow  can  also  be  significantly  modified  by  placing  a 
small  physical  obstruction  upstream  of  boundary  layer  separation  (e.g.,  Prandtl  et  al.,  1914).  The 
obstruction  (commonly  referred  to  as  a  “trip”)  is  thought  to  induce  premature  transition  to 
turbulence  within  the  boundary  layer,  and  consequently  delay  separation  presumably  owing  to 
enhanced  momentum  transfer  with  the  (outer)  cross  flow  (e.g.,  Fage  &  Warsap,  1929,  Igarashi, 
1986  and  Hover,  Tvedt  &  Triantafyllou,  2001).  In  fact,  similar  to  the  effects  of  a  synthetic  jet 
actuator,  the  asymmetric  placement  of  a  trip  wire  on  the  surface  of  a  circular  cylinder  can  lead  to 
reduction  in  drag  and  even  generate  lift  (James  &  Troung,  1972). 

The  present  report  describes  the  modifications  of  the  aerodynamic  characteristics  of  a  circular 
cylinder  at  sub-critical  Reynolds  numbers  using  integrated  synthetic  jet  actuators  and  miniature. 
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surface-mounted  passive  obstructions.  Even  though  the  details  of  the  interaction  dornains  of  the 
jet  and  the  obstruction  with  the  cross  flow  are  not  the  same,  it  is  shown  that  within  a  narrow 
range  of  operating  conditions,  a  properly  scaled  obstruction  can  induce  the  same  global 
aerodynamic  forces  as  an  active  fluidic  element.  It  appears  that  the  effect  of  these  Interaction 
domains  is  similar  to  that  of  circulation  bubbles  that  have  been  observed  upstream  of  separation 
during  the  (natural)  drag  crisis  at  the  critical  Reynolds  number  (e.g.,  Achenbach,  1968).  It  is 
remarkable  that  when  these  bubbles  are  formed  on  both  sides  of  the  cylinder  (near  ^=±100°), 
boundary  layer  separation  occurs  farther  downstream  at  (ftp«  147°)  and  the  drag  coefficient  is 
reduced  to  Cd  =  0.2  (Schewe,  1983).  These  changes  are  also  accompanied  by  a  substantially 
thinner  wake,  and  the  absence  of  vortex  shedding.  When  these  bubbles  ultimately  vanish  at 
supercritical  Reynolds  numbers  (when  the  boundary  layer  becomes  fully  turbulent),  the  drag 
coefficient  begins  to  increase  again  and  the  separation  point  moves  upstream. 

The  present  experiments  are  conducted  in  two  wind  tunnels  at  the  same  Reynolds  numbers  (for 
the  experiments  reported  here,  Reo  =  75,500).  Aerodynamic  measurements  are  conducted  in  an 
open-return  wind  tunnel  having  a  square  (0.91  m)  test  section,  while  high-resolution  PIV 
measurements  (8.88  pm/pixel,  and  spatial  resolution  of  0.133  mm)  along  tfie  surface  of  the 
cylinder  are  obtained  in  a  close-return  2-D  wind  tunnel  having  a  test  section  which  is  essentially 
a  5  cm  wide  slice  of  the  test  section  in  the  open  return  tunnel.  The  cylinder  model  (D  =  63.2 
mm)  can  be  rotated  about  its  spanwise  axis,  has  a  pair  of  adjacent  spanwise  jet  actuators  (2.5  inm 
apart,  each  0.5  mm  wide)  and  an  azimuthal  array  of  pressure  ports.  For  simplicity,  the  passive 
obstruction  is  a  surface-mounted  circular  cylinder  (diameter  d). 


1.3.1.  Global  Flow  Changes 

The  effects  induced  by  the  interaction  domain  between  the  jet  and  the  cross  flow  are  compared  to 
the  displacement  of  the  outer  flow  by  a  miniature  obstruction  having  a  circular  cross  section  that 
is  mounted  on  the  surface  of  the  cylinder  at  the  same  azimuthal  positions  of  the  jet.  The 
variations  of  Cl  and  for  various  obstruction  azimuthal  positions  and  diameters  are  shown  in 
Figures  I.3a  and  b,  respectively.  The  variation  of  the  forces  with  the  azimuthal  position  of  an 
obstruction  (d^D  =  0.0\7)  is  shown  in  Figure  I.3a.  Static  pressure  distributions  (not  shown) 
reveal  a  local  pressure  minimum  immediately  downstream  of  the  obstruction  (for  19°  <  /<  60°) 
that  is  associated  with  the  acceleration  of  cross  flow  fluid  over  the  separation  bubble  (cf.  Figure 
I.2b).  The  drag  is  reduced  by  a  significant  increase  in  the  base  pressure  {CJCp  reaches  a 
maximum  of  0.6  at  y=  50°)  and  then  abruptly  decreases  when  the  flow  separates.  As  shown  in 
Figure  I.3b,  the  characteristic  height  of  the  obstruction  (at  y=  60°)  has  a  profound  effect  on  the 
aerodynamic  performance  due  to  the  changes  in  the  size  of  the  separation  bubble.  There  is 
clearly  an  optimal  characteristic  scale  for  the  obstruction  for  which  the  lift  and  reduction  in 
pressure  drag  are  maximized.  As  the  size  of  the  obstruction  is  increased  beyond  d/D  >  0.02,  the 
separated  flow  does  not  reattach.  For  comparison.  Figure  I.3c  shows  the  variation  of  Cl  and 
Ci)=  AC[/Cdo  with  the  momentum  coefficient  of  a  synthetic  jet  that  is  placed  at  y=  60°  (the 
actuation  frequency  is  Sto  -  4.0).  While  both  the  passive  obstruction  and  the  jet  result  in  similar 


Figure  1.3  Variation  of  Cl  and  with:  a)  azimuthal  position  for  d/D  =  0.017,  b)  obstruction  diameter  at 
y=  60  °  and  c)  jet  momentum  coefficient  at  y=  60  ® 
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modifications  of  the  surface  pressure  distributions,  the  aerodynamic  forces  are  much  less 
sensitive  to  variations  in  momentum  coefficient  compared  to  changes  in  obstruction  size. 

1.3.2.  “Matched”  Aerodynamic  Performance  of  Jet  and  Passive  Obstruction 

The  details  of  the  interaction  domains 
of  the  passive  obstruction  and  the  jet 
are  compared  when  the  global 
aerodynamic  forces  are  approximately 
equal  (within  2%).  This  “matching”  is 
achieved  using  a  d/D  =  0.014 
obstruction  and  a  synthetic  jet  at  Sto  = 

4.0,  C^=  10’^  (Ci  =  0.49,  Cd  =  0.86 
for  the  obstruction  and  Cl  =  0.50, 

Cd  =  0.87  for  jet  actuation). 

Azimuthal  pressure  distributions  for 
the  “matched”  flows  (not  shown) 
reveal  that  aside  from  details  of  the 
pressure  peaks  (the  obstruction  yields 
a  narrower,  stronger  suction  peak 
compared  to  the  jet),  the  pressure  p® 
distributions  are  otherwise  very 
similar,  including  the  changes  in  base 
pressure  and  the  movement  of  the 
separation  points. 

The  interaction  between  the  actuation  Figure  1.4  Time-averaged,  normalized  spanwise  vorticity:  a) 
(jet  and  obstruction)  and  the  cross  baseline,  b)  obstruction,  and  c)  synthetic  jet  actuation. 

flow  on  the  top  surface  of  the  cylinder  Triangles:  closed-actuator  location,  open-separation  location. 

is  compared  with  the  baseline  flow 

using  maps  of  spanwise  vorticity  concentrations  that  are  obtained  from  PIV  measurements 
(Figures  I.4a-c).  The  baseline  flow  (Figure  I.4a)  separates  at  d~  95“  and  the  separated  flow 
domain  is  marked  by  the  presence  of  negative  (i.e.,  CCW)  vorticity  near  the  surface.  When  the 
obstruction  is  placed  on  the  surface 
(y=60“.  Figure  I.4b),  the  flow  separates 
locally  and  reattaches  at  0~15^  (i.e., 

-0.25  <  x/D  <  -0. 1 2)  although  the 
boundary  layer  downstream  of  the 
attachment  region  is  considerably  thicker 
than  for  the  baseline  flow  at  the  same 
azimuthal  domain.  Separation  occurs  at 
dsp=\\0°  and  the  streamwise  spreading 
rate  of  the  separated  shear  layer  appears  to 
be  smaller  than  for  the  unforced  flow 
presumably  owing  to  increased  speed  and 
turning  of  the  outer  flow.  The  vorticity 
maps  even  capture  the  blockage  effects 
and  thickening  of  the  boundary  layer 
upstream  of  the  obstruction.  Finally,  the 
effects  of  jet  actuation  are  shown  in  Figure 

I.4c.  It  is  remarkable  that  even  though  Figure  1.5  Cross  stream  velocity  (a)  and  Reynolds 

both  the  jet  and  the  obstruction  yield  stresses  (b).  a)  b-l  (o),  obst  (a),  jet  (•);  b)  u',u'Jul  - h-l 

similar  aerodynamic  forces  and  lead  to  (---),  obst  (•),  jet  (o):uyJui:  b-l  (—)  .  obst  (M),  jet 
flow  attachment  over  similar  azimuthal  (\3);';^fu^:b-l(—),obst(A),jet(A). 
extent  {dsp=\\\°  for  the  jet),  the  most  ' " 
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prominent  difference  between  the  vorticity  maps  is  that  the  flow  displacement  by  the  obstruction 
is  substantially  larger  than  the  displacement  by  the  interaction  domain  of  the  jet  (and  therefore 
the  suction  peak  induced  by  the  obstruction  is  larger).  Furthermore,  the  boundary  layer 
downstream  of  the  jet  interaction  domain  is  substantially  thinner  than  the  corresponding 
boundary  layer  downstream  of  the  obstruction.  Cross  stream  distributions  of  the  time-averaged 
azimuthal  velocity  and  Reynolds  stresses  downstream  from  the  obstruction  and  the  jet  (0=  100°) 
are  shown  in  Figures  3a  and  b,  respectively.  As  shown  in  Figure  I.5a,  at  this  azimuthal  position 
the  base  flow  is  already  separated  and  the  boundary  layer  downstream  of  the  obstruction  is 
significantly  thicker  and  the  velocity  outside  of  the  layer  is  somewhat  lower  than  in  the  boundary 
layer  induced  by  the  jet.  It  is  remarkable  that  the  distributions  of  Reynolds  stresses  downstream 
of  the  interaction  domains  (Figure  I.5b)  are  significantly  lower  than  in  the  baseline  flow  even 
when  the  flows  are  completely  attached  (not  shown).  That  this  reduction  is  also  prevalent 
(radially)  outside  of  the  boundary  layer,  indicates  that  the  two  interaction  domains  lead  to 
significant  reductions  in  velocity  fluctuations  that  are  associated  with  vortex  shedding  into  the 
wake.  In  a  similar  manner,  the  appearance  of  a  natural  separation  bubble  on  one  or  both  sides  of 
the  cylinder  at  the  critical  Reynolds  number  that  corresponds  to  the  drag  crisis  (Achenbach,  1968, 
Bearman,  1969)  results  in  substantial  changes  in  the  shedding  of  organized  vorticity  into  the 
wake.  These  findings  are  consistent  with  the  earlier  observations  of  Amitay,  Smith  &  Glezer 
(1998)  who  showed  that  the  Reynolds  stresses  in  the  wake  are  reduced  by  actuation  and  that  the 
spectral  peak  at  the  vortex  shedding  frequency  is  significantly  suppressed. 

1.4.  THE  INTERACTION  OF  A  SYNTHETIC  JET  WITH  A  “TRIPPED”  BOUNDARY  LAYER 

The  effects  of  synthetic  jet  actuation  in  a  cylinder  boundary  layer  that  is  presumably  “tripped”  by 
a  conventional  trip  wire  is  investigated  to  assess  whether  the  tripped  boundary  layer  impedes  the 
efficacy  of  the  jet  actuation.  The 
“tripped”  boundary  layer  is  turbulent-like 
flow  that  is  formed  by  symmetrically 
placing  a  pair  of  cylindrical  obstructions 
(d/D  =  0.02)  on  the  cylinder  surface  well 
at  /=  ±  34°  that  yield  a  turbulent-like 
(supercritical)  pressure  distribution  (as 
may  be  characterized  by  increased 
suction  peaks  and  base  pressure).  The 
changes  in  the  pressure  distribution  result 
in  a  substantial  reduction  (41%)  in 
pressure  drag  (Cd  =  0.82  compared  to  1 .4 
for  base  flow).  In  the  presence  of  the 
trip,  the  jet  is  more  effective  farther 
downstream  (near  }'=  100°,  Stp  =  4  and 
Cfi  =  3.S  X  10  ’).  The  asymmetry  induces 
additional  increase  and  decrease  in  the 
magnitude  of  the  suction  peak  on  the 
actuated  and  unactuated  surfaces  of  the 
cylinder  owing  to  the  change  in 
circulation  which  yields  Cl  =  0.36.  The 
increase  in  the  base  pressure  contributes 
to  a  further  reduction  in  drag  (Co  =  0.59) 
relative  to  the  “tripped”  flow.  Maps  of 
normalized,  time-averaged  spanwise 
vorticity  concentrations  with 
superimposed  cross  stream  velocity 
distributions  of  the  flow  fields  (i.e.,  the 
baseline,  and  in  the  absence  and  presence 
of  jet  actuation)  are  shown  in  Figures 


n/D 

Figure  L6  Normalized  spanwise  vorticity,  ReD  =  75,000: 
a)  baseline,  b)  tripped,  and  c)  actuated.  Boundary  layer 
separation  marked  with  open  triangles,  synthetic  jet 
actuator  location  marked  bv  closed  triangle. 
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I.6a-c.  Although  the  “tripped”  flow 
appears  to  be  attached  farther 
downstream  than  the  base  flow  and  the 
separating  shear  layer  appears  to  be 
deflected  towards  the  surface  of  the 
cylinder,  it  is  remarkable  that  boundary 
layer  separation  (as  determined  from 
zero  wall  shear  stress,  not  shown)  still 
occurs  at  nominally  6sp  =  93°  for  both 
the  baseline  and  tripped  flows. 

Corresponding  cross  stream  velocity 
distributions  at  ^  =  100°  (Figure  1.7) 
show  that  the  tripped  flow  is  indeed 
deflected  towards  the  wall  and  that  the  separated  flow  domain  occurs  within  a  narrow  cross 
stream  domain  that  is  less  than  one  millimeter  above  the  surface.  The  presence  of  separation 
bubbles  downstream  of  the  obstructions  at  y=  +  34°  apparently  modifies  the  local  streamwise 
pressure  gradient  and  results  in  a  thinner  boundary  layer  at  separation.  The  fact  that  the  location 
of  boundary  layer  separation  is  virtually  identical  to  that  of  the  baseline  flow  further  suggests  that 
transition  to  turbulence  is  incidental  in  the  alteration  of  the  global  aerodynamic  characteristics  of 
the  tripped  flow.  Nevertheless,  these  seemingly  minor  changes  in  the  flow  are  responsible  for  the 
substantial  changes  in  the  aerodynamic  forces  (incidentally,  a  similar  effect  is  observed  when  the 
synthetic  jet  is  positioned  at  y=  30°).  It  seems  plausible  that  a  similar  mechanism  is  associated 
with  the  onset  of  the  drag  crisis  at  the  critical  Reynolds  number  owing  to  the  appearance  of 
recirculating  flow  bubbles  near  ^  =  +  100°  (Achenbach,  1968).  Jet  actuation  of  the  “tripped” 
flow  causes  more  dramatic  changes  in  the  velocity  and  vorticity  distributions  which  include 
displacement  of  the  time-averaged  separation  point  to  dsp  =116°  and  a  downward  deflection  of 
the  separated  shear  layer.  These  changes  lead  to  a  thinner  wake  that  is  displaced  downward,  as 
noted  by  Glezer,  et  al.  (2003). 


Figure  1.7  Cross  stream  velocity  distributions,  6  =  100°: 
baseline  (O),  tripped  (0),  actuated  (A). 


Cross  stream  distributions 
of  the  Reynolds  stresses 
for  the  three  flows  are 
shown  in  Figures  I.8a-i.  In 
the  baseline  flow  (Figure 
I.8a-c),  the  three  stress 
components  increase  in 
magnitude  with 

downstream  distance,  and 
do  not  completely  vanish 
outside  of  the  boundary 
layer  even  upstream  of 
separation,  ostensibly  as  a 
result  of  potential 

fluctuations  that  are 
associated  with  the  global 
vortex  shedding.  As  noted 
above,  the  addition  of  the 
upstream  “trips”  results  in 
a  reduction  in  all  three 
components  of  Reynolds 
stress  across  the  boundary 
layer.  This  reduction, 

which  extends  outside  the 
layer  suggests  that  the  bulk 


Figure  L8  Normalized  Reynolds  Stresses:  (a,  d,  g),  UgUglU^  (b, 

e,  h),  u'lXlglul,  (e,f,  i);  baseline  (a,  b,  c),  tripped  (d,  e,J),  actuated  (g,  h,  i). 
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motion  due  to  global  vortex  shedding  is  diminished.  In  fact,  as  shown  in  Figures  I.8d-e  (tripped 
flow),  the  decrease  in  the  Reynolds  stress  persists  in  the  separated  shear  layer  (presumably  the 
result  of  reduction  in  both  global  potential  velocity  fluctuations  that  are  associated  with  vortex 
shedding  and  turbulence  production  within  the  separating  shear  layer  owing  to  stretching  in  the 
stain  field  of  the  large  scale  shed  vortices.  Jet  actuation  (Figures  I.8g-i)  brings  about  additional 
reduction  in  the  Reynolds  stresses  relative  to  the  “tripped”  flow,  both  within  the  boundary  layer 
and  the  separated  shear  layer  both  upstream  and  downstream  of  the  jet  actuator  (at  the  actuator 
orifice,  y=  100°,  the  peak  in  UgUg  has  approximately  the  same  magnitude  as  in  the  tripped  flow 
and  the  radial  velocity  fluctuations  are  larger  due  to  the  velocity  perturbations  introduced  by  the 
jet  itself.  However,  farther  downstream  the  turbulent  fluctuations  added  by  the  jet  are  rapidly 
dissipated  and  the  overall  levels  of  the  Reynolds  stresses  are  lower  than  in  the  “tripped”  flow 
suggesting  that  the  high  frequency  (small-scale)  actuation  that  are  introduced  by  the  jet,  lead  to 
enhanced  dissipation  in  the  base  flow  that  is  augmented  by  suppression  owing  to  the  presence  of 
global,  favorable  pressure  gradient. 

Even  though  it  is  difficult  to  assess  the  relative  importance  of  these  (potentially  competing) 
mechanisms  in  the  reduction  in  turbulence,  the  present  data  suggests  that  the  delay  in  separation 
and  the  concomitant  changes  in  the  global  aerodynamic  forces  are  not  merely  the  result  of 
turbulent  transition  within  the  surface  boundary  layer.  Rather,  these  changes  are  effected  by  the 
formation  of  the  interaction  domains  between  the  actuation  and  the  cross  flow  that  modify  the 
apparent  aerodynamic  shape  of  the  cylinder  (and  other  aerodynamic  surfaces). 
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II.  VIRTUAL  AERO-SHAPING  AT  LOW  ANGLES  OF  ATTACK 
11.1.  OVERVIEW 

The  pressure  drag  of  lifting  surfaces  at  low  angles  of  attack  is  reduced  with  minimal  lift  penalty 
by  fluidic  modification  of  the  apparent  aerodynamic  shape  of  the  surface.  The  pressure 
distribution  on  the  surface  is  modified  when  the  baseline  flow  is  fully  attached  by  the  formation 
of  a  small  stationary  recirculating  flow  domain  next  to  the  surface,  which  displaces  the  local 
streamlines  sufficiently  to  modify  the  local  pressure  distribution.  In  the  present  work,  the 
recirculating  flow  domain  is  formed  on  the  surface  of  a  clark-y  airfoil  by  combining  the 
activation  of  a  high-frequency  synthetic  jet  actuator  placed  downstream  from  a  miniature 
surface-mounted  passive  obstruction.  The  alteration  of  the  pressure  distribution  around  the 
airfoil  results  in  a  significant  reduction  in  pressure  drag  that  is  comparable  to  the  magnitude  of 
the  pressure  drag  of  the  baseline  configuration  with  minimal  loss  of  lift.  The  flow  structure  in 
the  vicinity  of  the  bubble  and  its  scaling  is  investigated  using  particle  image  velocimetry  (PIV). 

II.2  INTRODUCTION 

Much  of  the  recent  work  on  flow  control  techniques  with  the  objective  of  extending  the  post 
stall  flight  envelope  of  various  airfoil  configurations  has  focused  on  the  manipulation  of  flow 
separation  at  moderate  and  large  angles  of  attack  either  at  the  leading  edge  or  over  flaps  (e.g., 
Seifert  et  al.’).  This  has  been  typically  accomplished  by  exploiting  the  instability  of  the 
separating  shear  layer  and  its  receptivity  to  time-periodic  actuation  (e.g.,  pulsed  blowing)  on 
the  time  scale  of  the  flow  about  the  airfoil,  which  results  in  a  Coanda-like  unsteady 
reattachment.  Active  control  techniques  that  have  achieved  varying  degrees  of  separation 
control  by  manipulation  of  the  unstable  separated  free  shear  layer  have  included  external  and 
internal  acoustic  excitation  (e.g.,  Ahuja  and  Burrim,  and  Hsiao  et  al.  ),  vibrating  ribbons  or 
flaps  (e.g.,  Neuburger  and  Wygnanski'*),  and  steady  and  unsteady  blowing  or  bleed  (e.g., 
Williams  et  al.^,  and  Chang  et  al.^).  In  these  experiments,  the  time-periodic  actuation  was 
typically  applied  at  a  dimensionless  (reduced)  frequency,  F^~0(1)  such  that  the  actuation 
period  scaled  with  the  time  of  flight  over  the  length  of  the  reattached  flow. 

More  recently,  Smith  et  al.’  and  Amitay  et  al.®’  ^  demonstrated  the  suppression  of  separation 
over  an  unconventional  airfoil  at  moderate  Reynolds  numbers  (up  to  10  )  that  resulted  in  a 
dramatic  increase  in  lift  and  a  corresponding  decrease  in  pressure  drag.  Actuation  was  effected 
using  synthetic  (zero  mass  flux)  jet  actuators,  which  were  deliberately  operated  at  frequencies 
that  were  typically  an  order  of  magnitude  higher  than  the  characteristic  (shedding)  frequency  of 
the  airfoil  [i.e.,  ~  0(10)  rather  than  F"^  ~  0(1)].  These  authors  argued  that  the  interaction  of 

high-frequency  zero  net  mass  flux  jets  with  the  cross  flow  leads  to  local  modification  of  the 
apparent  aerodynamic  shape  of  the  flow  surface,  and,  as  a  result,  to  full  or  partial  suppression 
of  flow  separation.  Moreover,  the  recent  experiments  of  Erk‘°  demonstrated  the  suppression  of 
separation  on  an  FX61-184  airfoil  at  Reynolds  numbers  up  to  3x10  using  synthetic  jet 
actuation  at  frequencies  up  to  F"^  ~  0(100). 

The  interaction  between  a  cross  flow  over  a  2-D  circular  cylinder  and^  surface-mounted 
synthetic  jet  actuators  was  recently  investigated  in  detail  by  Honohan  et  al.  .  These  authors 
showed  that  when  the  jets  are  operated  on  a  time  scale  that  is  well  below  the  characteristic  time 
scale  of  the  base  flow,  their  interaction  with  the  cross  flow  leads  to  the  formation  of  distinct 
quasi-steady  flow  regions  near  the  surface  and  displacement  of  the  local  streamlines  induces  an 
apparent  or  virtual  change  in  the  shape  of  the  surface  and  in  the  local  pressure  gradient.  Tlie 
acceleration  of  the  cross  flow  around  the  interaction  domain  is  accompanied  by  substantial 
alterations  of  the  streamwise  pressure  gradient  both  locally  and  globally.  As  a  result  the 
surface  boundary  layer  downstream  of  the  interaction  domain  becomes  thinner  allowing  the 
flow  to  overcome  stronger  adverse  pressure  gradients  and  therefore  delaying  (or  altogether 
suppressing)  flow  separation. 
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The  concept  of  modifying  the  apparent  aerodynamic  shape  of  aero-surfaces  in  order  to 
prescribe  the  streamwise  pressure  distribution  and  therefore  to  influence  its  aerodynamic 
performance  is  not  new  and  was  addressed  in  a  substantial  body  of  work  in  the  40s  and  50s. 
One  example  of  this  approach  is  the  use  of  a  stationary,  trapped  vortex  to  alter  the  apparent 
local  surface  curvature  and  therefore  the  direction  of  the  external  flow.  Perkins  and  Hazen 
used  a  suction  slot  near  the  trailing  edge  of  an  airfoil  and  an  upstream  cavity  to  trap  a  vortex 
near  the  trailing  edge  and  reported  substantial  lift  (Cl  =  1 .2)  at  zero  angle  of  attack  at  relatively 
low  (3%)  suction  coefficient.  Mandl'^  used  a  split  flap  located  on  the  pressure  side  of  an  airfoil 
combined  with  suction  (at  the  flap)  to  trap  a  vortex  and  to  enforce  tangential  flow  near  the 
trailing  edges  of  the  airfoil  and  the  flap  and  alter  the  Kutta  condition  of  the  airfoil. 

The  present  work  focuses  on  the  modification  of  the  apparent  aerodynamic  shape  of  an  airfoil 
at  low  angles  of  attack  when  the  baseline  flow  is  fully  attached  with  the  objective  of  reducing 
the  pressure  drag  with  minimal  or  no  penalty  in  lift  Chatlynne  et  al.’^  and  Amitay  et  al.  ).  This 
is  achieved  by  the  formation  of  a  stationaiy  recirculating  flow  domain  next  to  the  surface  that 
alters  the  flow  above  the  airfoil  by  inducing  a  displacement  of  the  local  streamlines.  While 
such  a  recirculating  domain  can  be  formed  in  principle  by  an  array  of  Jet  actuators,  in  the 
present  work  this  domain  is  formed  by  combining  the  activation  of  a  high-frequency  synthetic 
jet  actuator  that  is  placed  downstream  from  a  miniature  surface-mounted  passive  obstruction. 

11.3  EXPERIMENTAL  APPARATUS  AND  PROCEDURE 

Fluidic  modification  of  the  aerodynamic  performance  of  lifting  surfaces  is  investigated  using  a 
2-D  Clark-Y  airfoil  having  a  chord  measuring  229  mm  and  thickness  to  chord  ratio  t/c  of  0.24. 
Control  is  effected  using  an  array  of  four  collinear  spanwise  rectangular  synthetic  jet  actuators 
that  are  integrated  into  the  center  section  of  the  airfoil  as  shown  schematically  in  Figure  II.4. 
Each  actuator  jet  has  a  flush-mounted  140  mm  long  and  0.5  mm  wide  (in  the  streamwise 
direction)  orifice.  The  jets  issue  normal  to  the  surface  of  the  airfoil  and  are  located  at 
x/c  =  0.22,  0.28,  0.36,  and  0.44  (measured  from  the  leading  edge).  Streamwise  fences  are 
placed  near  the  (spanwise)  edges  of  the  actuator  orifices  to  maintain  a  nominally  two- 
dimensional  flow  in  the  controlled  section  of  the  airfoil.  The  actuator  jets  are  driven  by 
piezoelectric  membranes  built  into  a  cavity  underneath  the  skin  of  the  airfoil.  The  airfoil  is  also 
instrumented  with  a  circumferential  array  of  41  pressure  ports  (located  at  mid-span)  that  are 
each  connected  to  an  external  high-speed  pressure  measurement  system.  The  airfoil  is  mounted 
in  a  low-speed  open-return  wind  tunnel  having  a  square  test  section  measuring  91  cm  on  a  side 
and  maximum  air  speed  of  32  m/s  (the  free-stream  turbulence  level  is  lower  than  0.25%).  The 
upper  and  lower  walls  of  the  wind  tunnel  test  section  are  adjusted  to  compensate  for  blockage 
created  by  the  airfoil. 

In  the  present  work,  the  performance  of  the  synthetic  jet  actuators  is  measured  using  their 
momentum  coefficient 


C„  =- 


I, 


\PoVlc 


where  7 j  is  the  time-averaged  momentum  flux  per  unit  length  during  the  outstroke  and  is  given 
by 

r 

Jj  =^Pjb\<u]{(p)>d^, 


X  =  T/2  (T  is  the  period  of  the  diaphragm  motion),  pj  and  po  are  the  jet  and  free-stream  fluid 
densities,  respectively,  b  is  the  jet  orifice  width,  c  is  the  chord,  Uo  is  the  free-stream  velocity, 
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and  <Uj((|))>  is  the  phase-averaged  velocity  at  the  jet  exit  plane.  During  calibration,  the  orifice 
velocity  of  the  jets  is  measured  using  a  miniature  hot-wire  sensor  following  the  procedure  of 
Smith  and  Glezer’'*. 

The  velocity  and  vorticity  fields  in  the  cross-stream  (x-y)  plane,  z  =  0,  above  the  airfoil  (i.e.,  on 
the  suction  side)  are  measured  using  particle  image  velocimetry  (PIV).  The  flow  is  seeded  with 
sub  micron  smoke  particles  and  it  is  illuminated  using  a  double-pulse  ND-YAG  laser.  Image 
pairs  are  captured  using  a  1008x1016  pixel  CCD  camera  with  a  magnification  of  21  and 
86|im/pixel  for  global  and  detailed  views,  respectively.  In  each  image,  the  velocity  vectors  are 
computed  using  a  standard  cross-correlation  technique.  A  given  PIV  map  is  averaged  over  300 
realizations  (image  pairs). 

11.4.  AERODYNAMIC  MODIFICATION  OF  THE  AIRFOIL  PERFORMANCE 

Although  the  primary  objective  of  the  present  work  is  to  demonstrate  the  effect  of  fiuidic 
actuation  on  the  aerodynamic  performance  of  the  airfoil  at  low  angles  of  attack  (i.e.,  in  the 
absence  of  stall),  it  is  instructive  to  gauge  the  effectiveness  of  the  present  actuators  on  the 
baseline  configuration  over  a  broad  range  of  angles  of  attack.  To  this  end,  the  pressure 
distribution  about  the  airfoil  was  measured  for  a  between  -1  and  24°  in  the  presence  and 
absence  of  actuation.  Figure  Il.la-c  show  the  distributions  of  the  pressure  coefficient  for 
a=  1.5°,  9°  and  19°,  respectively.  The  chord  Reynolds  number  is  Rec  =  250,000,  the  active 
actuator  jet  is  located  at  x/c  =  0.22  and  the  jet  momentum  coefficient  is  C^  =  1.9 10'^.  As 
evident  from  the  data,  at  a  =  1.5®  (Figure  II.  la),  the  baseline  flow  in  the  absence  of  the  control 
jet  is  fully  attached  over  the  entire  surface  of  the  airfoil.  It  is  noteworthy  that  when  the  jet  is 
activated,  the  suction  peak  decreases  somewhat  (by  2%)  and  moves  slightly  downstream  and 
the  (suction)  magnitude  of  the  pressure  coefficient  remains  lower  than  in  the  absence  of 
actuation  through  the  trailing  edge  of  the  airfoil.  These  changes  in  Cp  lead  to  a  small  decrease 
in  lift  (8%,  as  shown  in  Figure  II.2a  above)  while  the  pressure  drag  coefficient  remains 
unchanged.  This  translates  to  a  reduction  in  the  lift  to  pressure)  drag  ratio.  The  present 
measurements  show  that  in  the  absence  of  control,  the  baseline  airfoil  stalls  at  oc «  7°  (cf.  Figure 
II.2b).  In  Figure  II.  lb,  (a  =  9°)  the  baseline  flow  is  stalled  and  when  the  jet  (which  is  located  in 
the  separated  flow  domain)  is  activated,  the  flow  becomes  fully  attached  along  most  of  the 
surface  of  the  airfoil  (except  possibly  the  streamwise  domain  x/c  >  0.75)  and  the  resulting 
pressure  distribution  exhibits  a  large  suction  peak  (Cp  =  -2.8)  near  x/c  =  0.2,  which  is  near  the 
point  of  maximum  thickness  (x/c  =  0.25).  The  apparent  increase  in  circulation  about  the 
attached  airfoil  affects  the  pressure  distribution  around  the  entire  airfoil  and  in  particular  it 
appears  that  the  front  stagnation  point  moves  towards  the  pressure  side  and  that  the  pressure 
increases  along  the  bottom  surface  of  the  airfoil.  The  effectiveness  of  the  actuation  diminishes 


'x/c  '  x/c  x/c 

Figure  ILL  Pressure  distribution  around  the  aiifoil  at  a  =  1.5"  (a),  9®  (b)  and  19^  (c).  Re^  =  250,000  and 
Cf,  =  1.9x1  (  o  )  baseline  and  (  •  )  actuation. 
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at  high  angles  of  attack.  At  a  =  19°  (Figure  11.  Ic),  the  pressure  peak  is  lower  (Cp  =  -2.1)  and 
moves  farther  upstream  (x/c  =  0.1)  and  it  appears  that  the  flow  is  attached  only  through 
x/c  =  0.5  and  is  separated  farther  downstream.  Note  that  even  though  the  flow  is  partially 
separated,  the  pressure  distribution  along  the  pressure  side  is  almost  identical  to  the  distribution 
in  the  absence  of  actuation. 


Figure  11.2.  Lift  coefficient  (a),  pressure  drag  coefficient  (b)  and  lift  to  pressure  drag  ratio  (c)  as  a  function  of  a 
(  o  )  baseline  and  (  •  )  actuation. 


The  pressure  distributions  on  the  airfoil  in  the  presence  and  absence  of  actuation  are  integrated 
to  yield  the  lift  coefficient,  the  pressure  drag  coefficient  and  the  lift-to-pressure  drag  ratio 
(Figure  II.2a-c,  respectively).  As  noted  above,  these  data  show  that  in  the  absence  of  actuation 
the  airfoil  stalls  at  a  =  7°  and  thereafter  the  lift  remains  almost  invariant  with  angle  of  attack 


Figure  II.3.  Pitching  moment  coefficient  (a)  and  the 
center  of  pressure  (b)  as  a  junction  of  a  (  o  ) 
baseline  and  (  •  )  actuation. 


although  the  pressure  drag  continues  to 
increase.  The  stall  results  in  a  sharp  decrease 
(an  order  of  magnitude)  in  the  lift  to  pressure) 
drag  ratio.  As  noted  in  connection  with  Figure 
II. la  above,  actuation  at  low  angles  of  attack 
(i.e.,  pre-stall,  a  <  5°),  results  in  a  nominal  8% 
reduction  in  lift  while  the  drag  coefficient 
remains  the  same  and  therefore  the  net 
reduction  in  the  peak  level  of  L/Dp  (at  a  =  1 .5°) 
is  approximately  10%.  At  post  stall  angles  of 
attack,  the  actuation  leads  to  a  substantial 
augmentation  in  lift  that  is  accompanied  by  a 
reduction  in  the  (pressure)  drag  coefficient  and 
consequently  to  a  higher  lift-to-pressure  drag 
ratio.  In  the  presence  of  actuation,  the  length  of 
the  separated  flow  domain  upstream  of  the 
trailing  edge  increases  with  angle  of  attack  and 
there  is  a  relatively  small  decrease  in  the  peak 
suction  pressure  (e.g..  Figure  Il.lb  and  c).  The 
net  effect  is  that  in  the  post-stall  domain  for  the 
baseline  airfoil  in  the  absence  of  control,  the 
actuation-augmented  lift  remains  almost 
invariant  with  angle  of  attack  and  therefore  the 
gradual  increase  in  pressure  drag  results  in  a 
similar  decrease  in  L/Dp. 

The  effects  of  the  actuation  on  the  location  of 
the  airfoil’s  center  of  pressure  and  the  pitching 
moment  (computed  relative  to  c/4)  are  shown 
in  Figure  II.3a  and  b.  Stall  results  in  a 
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substantial  change  in  the  center  of  pressure  of  the  baseline  airfoil  (nominally  35%  at  a  =  9°), 
which  is  accompanied  by  a  (mild)  decrease  in  the  magnitude  of  the  pitching  moment.  The 
actuation  eliminates  the  abrupt  change  in  the  center  of  pressure  and  within  the  post  stall  domain 
shown  here,  Xcp  varies  by  no  more  than  10%  while  the  pitching  moment  decreases  compared  to 
the  baseline  configuration. 

•  •  -3  • 

The  data  shown  above,  suggests  that  relatively  low-level  actuation  [i.e.,  Cp  of  0(10' )]  is  most 
effective  at  post  stall  angles  of  attack  of  the  baseline  airfoil  but  that  the  induced  aerodynamic 
changes  at  low  angles  of  attack  are,  in  effect,  marginal.  Other  work  at  Georgia  Tech  (not 
shown)  has  demonstrated  that  control  effectiveness  can  be  realized  at  low  angles  of  attack  with 
increased  actuation  levels  (or  higher  C^).  However,  as  indicated  by  the  present  results,  a  new 
approach  to  the  modification  of  the  aerodynamic  performance  of  airfoils  at  low  angles  of  attack 
may  be  necessary  if  it  is  also  desirable  to  maintain  low  actuation  levels. 

In  what  follows,  a  new  flow  control  technique 
is  presented  with  the  objective  of  reducing  the 
pressure  drag  at  low  angles  of  attack  with 
minimal  or  no  penalty  in  lift.  The  approach  is 
to  modify  the  pressure  distribution  on  the 
surface  of  the  airfoil  by  exploiting  virtual 
aero-shaping  of  the  surface  that  is  achieved  by 
the  formation  of  a  stationary  recirculating 
flow  domain  next  to  the  surface.  The 
formation  of  this  domain  alters  the  flow  above 
the  airfoil  by  inducing  a  displacement  of  the 
local  streamlines  that  is  sufficient  to  modify 
the  pressure  distribution  on  the  surface  of  the 
airfoil  even  at  small  angles  of  attack.  While  such  a  recirculating  domain  can  be  formed  in 
principle  by  an  array  of  jet  actuators,  in  the  present  work  this  domain  is  formed  by  combining 
the  activation  of  a  high-frequency  synthetic  jet  actuator  that  is  placed  downstream  from  a 
miniature  surface-mounted  passive  obstruction  (as  shown  schematically  in  Figure  II.4). 

In  a  proof  of  concept  experiment,  the 
characteristic  height  of  the  obstruction  (above 
the  surface  of  the  airfoil)  is  h  =  9.310'^c  (i.e., 
less  than  1%  of  the  chord)  and  it  is  placed  at 
Xo/c  =  0.19.  The  synthetic  jet  is  located  at 
Xj/c  =  0.22  (i.e.,  downstream  of  the 

obstruction)  and  is  operated  at  0^=  1.9 10'^. 

In  these  experiments,  the  angle  of  attack  is  1° 
and  Rec  =  381,000.  The  distribution  of  the 
pressure  coefficient  around  the  circumference 
of  the  airfoil  for  the  baseline  configuration 
and  in  the  presence  of  actuation  are  shown  in 
Figure  II.5.  At  this  angle  of  attack  the 
baseline  flow  is  attached  to  the  surface  of  the 
airfoil.  The  activation  of  the  synthetic  jet 
actuator  downstream  of  the  obstruction  results 
in  an  increase  (42%)  in  the  magnitude  of  the 
suction  peak  and  its  displacement  in  the 
upstream  direction  (to  x/c  =  0.21).  The  pressure  recovery  is  such  that  magnitude  of  the  static 
pressure  for  x/c  >  0.33  is  smaller  than  for  the  baseline  configuration  (i.e.,  within  this  domain, 
the  pressure  difference  between  the  suction  and  pressure  sides  is  smaller  when  the  control  is 
applied)  indicating  that  there  may  be  a  net  reduction  in  pressure  drag.  Note  that  upstream  of 
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Figure  IL5.  Pressure  distribution  around  the  airfoil 
ata=I.  Rec  =  381,000  and  C^  =  I.9xl(r\  (  o  ) 
baseline  and  (  m  )  actuation. 


the  obstruction  the  pressure  is  slightly  lower  than  for  the  baseline  configuration  (i.e.,  less 
negative)  due  to  blockage  caused  by  the  obstruction. 

These  effects  are  evident  from  distributions  of  the  lift,  pressure  drag,  and  the  pressure  drag 
polar  (Figures  6a-c)  that  are  computed  from  a  series  of  pressure  measurements  over  a  range  of 


Figure  II.6.  Lift  coefficient  (a)  and  pressure  drag  coefficient  (b)  as  a  function  of  ce,  and  the  pressure  drag  polar  (c) 
(  o  )  baseline  and  (  •  )  actuation. 


0.3 


y/c 
0.2 

Figure  IL  7.  Time-averaged  velocity  vector  maps  above  the  airfoil  Global  views:  baseline  (a)  and  actuation  (b), 
zoom  in:  baseline  (c)  and  actuation  (d). 
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angles  of  attack  (-9®<a<5°).  For  this  range  of  angles  of  attack,  activation  of  the  control 
results  in  reductions  in  both  of  (pressure)  drag  and  the  lift  coefficients  (the  reduction  in  lift  is 
almost  uniform  throughout  the  range  of  measurements).  However,  while  the  reduction  in  the 
lift  is  about  10%  with  respect  to  the  baseline  confi^ration,  the  corresponding  reduction  in  the 
(pressure)  drag  is  greater  than  30%  resulting  in  an  increase  of  the  lift-to-pressure  drag  ratio  for 
a  given  angle  of  attack.  Moreover,  the  pressure  drag  polar  (Figure  6c)  shows  that  for  Cl  <  0.25, 
the  pressure  drag  of  the  modified  airfoil  decreases  significantly.  The  reduction  in  (pressure) 
drag  without  a  substantial  loss  in  lift  may  be  important  in  cruise  conditions  and  lead  to 
significant  savings  in  fuel  over  long  distances. 

The  effect  of  the  actuation  on  the  flow  field  above  the  surface  of  the  airfoil  is  investigated  using 
particle  image  velocimetry  (PIV)  in  the  x-y  plane  (z  =  0).  Each  image  is  comprised  of  three 
partially  overlapping  frames  measuring  10cm  x  10cm  and  consists  of  300  realization  pairs. 
Figures  II.7a-d  show  cross-stream  maps  of  the  time-averaged  velocity  vector  of  the  baseline 
configuration  (Figures  II. 7a  and  c)  and  of  the  controlled  flow  (Figures  II.7b  and  d).  Magnified 
views  of  the  domain  of  interaction  between  the  control  actuation  and  the  cross  flow  are  shown 
in  Figures  II.7c  and  d  (the  magnification  in  Figures  II.7c  and  d  is  21}im/pixel  compared  to 
86|xm/pixel  in  Figures  II.7a  and  b).  In  the  absence  of  control,  the  time-averaged  velocity  field 
is  clearly  attached  to  the  surface  of  the  airfoil  (Figures  II.7a  and  c).  When  actuation  is  applied 
(Figure  II.7b),  the  velocity  vector  field  above  the  surface  of  the  airfoil  appears  to  be 
qualitatively  similar  to  the  corresponding  velocity  field  in  the  absence  of  the  control  (Figure 
II.7a).  Although  the  surface  boundary  layer  is  thicker  near  the  trailing  edge  of  the  airfoil  in  the 
controlled  flow,  ostensibly  as  a  result  of  a  sub  optimal  pressure  recovery.  The  high- 
magnification  (zoomed  in)  images  of  the  flow  in  the  vicinity  of  the  actuation  jet  show  the 
formation  of  a  closed  recirculation  bubble  downstream  of  the  obstruction  (the  location  of  the 
active  jet  is  marked  with  an  arrow).  It  is  remarkable  that  despite  its  relatively  small 
characteristic  scales  (cross-stream  and  streamwise  dimensions  of  0.03c  and  0.1c,  respectively), 
the  presence  of  this  bubble  is  sufficient  to  effect  the  aerodynamic  modifications  depicted  in 
Figure  II.b.The  internal  structure  of  the  recirculation  bubble  is  investigated  using  phase- 
averaged  velocity  vector  maps.  The  vector  maps  in  Figures  II.8a-f  are  acquired  at  six  time 
increments  during  the  jet  actuation  cycle  where  Figures  II.8a  and  d  correspond  to  the  beginning 
of  the  blowing  and  suction  cycles.  At  the  beginning  of  the  blowing  cycle  (Figure  II.8a),  the 
recirculation  domain  includes  a  clockwise  (CW)  vortex  centered  approximately  at  x'  =  6.rh 
downstream  from  the  obstruction.  The  image  also  shows  the  blockage  upstream  of  the 
obstruction  and  the  separated  flow  downstream  from  its  cross-stream  edge.  A  smaller  CW 
vortex  is  visible  farther  (12.2  h)  downstream.  As  the  blowing  proceeds  (Figures  II.8b  and  c), 
the  upstream  (larger)  CW  vortex  seems  to  intensify  as  is  evidenced  by  the  relatively  strong 
downwash  that  is  induced  near  the  downstream  edge  of  the  primary  bubble  and  is  advected 
farther  downstream  (the  maximum  streamwise  length  of  the  primary  bubble  is  ~3.8  h).  These 
data  also  suggest  that  a  weak  CCW  vortex  is  formed  on  the  left  side  of  the  jet  orifice  (the 
nominally  2-D  jet  forms  a  counter  rotating  vortex  pair  during  the  blowing  cycle).  At  the  same 
time  it  appears  that  the  strength  of  the  secondary  (downstream)  vortex  is  somewhat  reduced. 
The  effect  of  the  suction  stroke  is  apparently  to  ingest  the  weak  CCW  vortex  (immediately 
downstream  of  the  obstruction)  and  to  induce  the  formation  of  a  second  CW  vortex  within  the 
primary  bubble  (Figure  II.8e),  which  intensifies  as  the  downstream  (earlier)  CW  vortex 
becomes  weaker.  By  the  time  the  cycle  repeats  Itself,  the  latter  vortex  is  advected  outside  of 
the  downstream  edge  of  the  primary  bubble.  This  sequence  of  images  clearly  indicates  that  the 
streamwise  extent  of  the  bubble  is  of  the  order  of  one  to  two  wavelengths  of  the  actuation 
frequency  (e.g.,  based  on  half  the  local  free  stream  velocity).  More  importantly,  similar  to  the 
earlier  observations  of  Honohan  et  al."  farther  downstream,  the  actuation  frequency  has 
negligible  temporal  effect  on  the  vorticity  distribution  within  the  wall  boundary  layer. 

The  characteristic  streamwise  scale  of  the  recirculation  bubble  (and  consequently  the 
aerodynamic  performance  of  the  airfoil)  can  be  extended  by  moving  the  actuation  jet  farther 
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Figure  II.8.  Phase-averaged  velocity  vector  maps  above  the  airfoil,  actuation  on.  t/T  =  0  (a),  0.29  (b),  0.42 
(c),  0.5  (d),  0.71  (e)  and  0.93  (f). 

downstream  relative  to  the  obstruction.  For  an  obstruction  having  a  characteristic  height 
(h/c  =  0.01)  that  is  placed  at  x/c  =  0.22,  the  jet  streamwise  position  (relative  to  the  obstruction) 
is  increased  from  x’/h  =  0  to  6  (x’  is  the  streamwise  distance  measured  from  the  obstruction). 
The  resulting  distributions  of  the  pressure  coefficient  are  shown  in  Figure  II.9  (the  distribution 
for  the  baseline  configuration  in  the  absence  of  actuation  is  also  shown  for  reference).  The 
most  prominent  feature  of  these  distributions  is  the  increase  in  the  magnitude  of  the  low 
(suction)  pressure  within  the  interaction  domain  of  the  actuation  with  the  cross  flow.  The 
upstream  jet  (x’/h  =  0)  results  in  a  large  suction  peak  at  x/c  =  0.25  (-2.43  compared  to  -1.55 
for  the  baseline  configuration)  followed  by  a  rapid  pressure  recovery  for  0.25  <  x/c  <  0.32,  and 
thereafter  the  pressure  recovery  towards  the  trailing  edge  is  more  gradual.  Compared  to  the 
upstream  jet,  activation  of  the  jet  that  is  located  at  (x’/h  =  6)  yields  a  lower  suction  peak 
(around  x/c  =  0.25)  and  a  wider  pressure  distribution  around  it.  However,  the  pressure  recover 
for  x/c  >  0.4  is  similar  for  both  cases  (note  that  the  pressure  in  the  presence  of  actuation  is 
lower  than  for  the  baseline  configuration  indicating  a  reduction  in  pressure  drag). 

The  corresponding  distributions  of  lift,  pressure  drag,  and  the  pressure  drag  polar  are  shown  in 
Figures  Il.lOa-c,  respectively.  When  the  control  jet  is  just  downstream  of  the  obstruction 
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(x’/h  =  0)  the  pressure  drag  coefficient  is  reduced  throughout  the  entire  measurement  range 
while  the  lift  is  almost  unchanged 


x/c 

Figure  11.9.  Pressure  distribution  around  the  aitf oil  at  a  =  I. ).  (  o  )  baseline,  (  •  )  upstream  jet  is  on  and 

(  *  )  downstream  Jet  is  on. 

compared  to  the  baseline  configuration.  In  fact,  the  pressure  drag  polar  shows  that  for  a  given 
magnitude  of  the  lift  the  (pressure)  drag  is  reduced  (the  major  effect  is  for  a<  1”).  The 
performance  of  the  actuation  (in  terms  of  reduction  in  pressure  drag  with  a  minimum  penalty  in 
lift)  is  somewhat  degraded  when  the  jet  is  placed  farther  downstream  (x’/h  =  6)  and  seems  to  be 
related  to  the  characteristic  dimensions  of  the  interaction  bubble.  It  is  noted  that  the 
downstream  jet  can  not  maintain  an  attached  flow  (in  the  presence  of  the  obstruction)  at  angles 
of  attack  above  3”. 

As  noted  above,  based  on  the  data  in  Figure  II.9  it  appears  that  placement  of  the  actuation  jet 
farther  downstream  extends  the  streamwise  length  of  the  recirculating  bubble  and  hence 
diminishes  the  magnitude  of  the  suction  peak  as  the  external  flow  turns  around  the  bubble. 
This  effect  is  visible  in  time-averaged  vector  maps  obtained  from  PIV  images  in  the  streamwise 
domain  0.15  <  x/c  <  0.4  (Figures  II.l  la  and  b).  These  data  show  that  the  maximum  height  of 
the  bubble  increases  from  1.07h  to  1.7  Ih  and  its  length  increases  from  6.5h  (0.066c)  for 
x’/h  =  0  to  12.9h  (0.1 6h).  These  data  confirm  that  the  pressure  distribution  in  the  vicinity  of 
the  bubble  can  be  adjusted  (and  controlled)  by  extending  its  streamwise  length  and  the  rate  of 
turn  of  the  outer  flow. 


Figure  ILIO.  Lift  coefficient  (a)  and  pressure  drag  coefficient  (b)  as  a  function  of  a,  and  the  pressure  drag 
polar  (c).  (  o  )  baseline,  (  •  )  upstream  jet  is  on  and  (  *  )  downstream  jet  is  on. 


In  the  experiments  discussed  above,  the  momentum  coefficient  of  the  jet  was  invariant 
(Cn  =  1 .9x10'^).  The  sensitivity  of  the  actuation  effectiveness  (as  measured  by  the  aerodynamic 
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forces)  to  is  measured  over  the  range  1.8 10'^  <  ^  1.9 10'  for  at  a=  1°,  X|/c  =  0.22.  The 

resulting  pressure  distributions  around  the  airfoil  (Figure  11.12)  are  very  similar  and  the 
differences  in  the  vicinity  of  the  actuation  domain  are  rather  small.  For  example,  an  increase  in 
Cn  from  3.110'^  to  1.9  10'^  results  in  an  increase  of  the  suction  peak  by  10%.  Therefore,  it 
appears  that  for  this  range  of  momentum  coefficients  the  size  of  the  bubble  is  almost  invariant. 
These  data  indicate  that  at  least  for  the  present  actuation  the  effect  of  on  the  lift  and 
(pressure)  drag  (not  shown)  and  thus  on  the  aerodynamic  performance  is  minimal. 


Figure  11.11.  Time-averaged  velocity  vector  maps  above  the  airfoil.  Upstream  jet  on  (a)  and  downstream  Jet  on  (b). 

Finally,  given  the  above  results,  it  might  be  argued  that  similar  effects  may  be  obtained  by 
simply  replacing  the  recirculating  flow  domain  that  results  from  the  actuation  with  a  solid 
extension  of  the  surface  having  the  same  shape.  This  notion  is  tested  by  comparing  the  effects 
of  the  actuation  and  of  a  passive  surface  modification  on  the  pressure  distribution  around  the 


Figure  11. 12.  Pressure  distribution  around  the  airfoil  at  a  = 
(  •  )C^  =  l.9xl(y\(  A  )5.4xl(r\(*  )3.lxia^and( 


f.  Global  distribution  (a)  and  detailed  distribution  (b). 
♦  )l.8.1xl(f\ 
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Figure  11.13.  Pressure  distribution  around  the  airfoil  at  a  =  1".  (  o  )  baseline,  (  •  )  actuation  and  ( - ) 

solid  surface  extension. 

airfoil  (Figure  13).  The  passive  extension  to  the  surface  was  fabricated  (using  stereo¬ 
lithography)  based  on  the  streamline  maps  that  were  computed  from  the  PIV  data.  This 
extension  extends  approximately  0.01c  above  the  surface  and  0.12c  in  the  streamwise  direction. 

As  shown  in  Figure  13,  the  surface  extension  leads  to  a  complete  flow  separation  ostensibly 
because  the  flow  above  the  airfoil  cannot  overcome  the  local  adverse  pressure  gradient  that  is 
imposed  by  the  solid  extension.  Furthermore,  it  is  clear  that  the  solid  extension  does  not 
duplicate  the  boundaiy  conditions  on  the  surface  of  the  recirculating  bubble  in  terms  of  slip  and 
penetration. 
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11.5  SCALING  OF  THE  INTERACTION  DOMAIN 

As  shown  above,  the  pressure  drag  at  low  angles  of  attack  can  be  reduced  with  minimal  penalty 
in  lift  by  forming  a  stationary,  recirculating  flow  domain  next  to  the  surface  and  thereby 
modifying  its  apparent  aerodynamic  shape  and  altering  the  surface  pressure  distribution  even  at 
small  angles  of  attack.  The  recirculating  domain  is  formed  by  combining  the  activation  of  a 
high-frequency  synthetic  jet  actuator  that  is  placed  downstream  from  a  miniature  surface- 
mounted  passive  obstruction  (as  shown  schematically  in  Figure  11.4). 

In  this  section,  attention  is  restricted  to  the  characteristic  streamwise  scale  of  the  recirculation 
domain  at  a  fixed  angle  of  attack  (a  =  3°).  This  domain  is  altered  by  varying  the  actuation 
frequency  and/or  the  jet  momentum  coefficient.  The  resulting  distributions  of  the  pressure 
coefficient  are  shown  in  Figures  II.14a-d  (the  distribution  for  the  baseline  configuration  in  the 
absence  of  actuation  is  also  shown  for  reference)  fory^c/=  139,  187,  362,  and  850  Hz, 
respectively.  The  most  striking  feature  of  these  pressure  distribution  is  the  appearance  of  a 
local  minimum  in  the  suction  pressure  just  upstream  of  the  obstruction  (e.g..  Figure  II.  14a  atx/c 
=  0.2).  This  minimum  vanishes  as  the  actuation  frequency  is  increases  and  is  virtually 
undetectable  at  850  Hz.  The  presence  of  this  minimum  suggests  a  blockage  effect  just 
upstream  of  the  obstruction  that  can  be  attributed  to  the  cross-stream  height  of  the  recirculating 
domain.  This  blockage  leads  to  a  local  reduction  in  suction  pressure  compared  to  the  baseline 
flow  and  therefore  to  loss  of  lift  and  potentially  an  increase  in  pressure  drag. 


Figure  11.14.  Pressure  distributions  around  the  airfoil  at  a=  3^  and  Rec  -  381,000  for  (a)  faa  -  I39Hz, 
(b)  187  Hz.  (c)  362  Hz  and  (d)  850  Hz.  (-^  baseline  and  actuated:  (•)  =  1.2x10-^ ,  (A>  3.  Zx/tT' ,  and 

(*)1.7xl(r\ 
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The  effect  of  the  blockage  on  the  flow  field  can  be  assessed  from  cross-stream  velocity  profiles 
that  are  measured  above  the  obstruction  (using  particle  image  velocimetry).  Figure  11.15  shows 
two  such  profiles  measured  at  x/c  =  0.16  and  0.21  (above  the  obstruction)  for_^c/  =139  HZ  and 
850  Hz  (the  jet  location  is  marked  with  an  arrow).  The  velocity  vectors  that  correspond  to  the 
high  frequency  actuation  not  only  exhibit  higher  speed  near  the  surface,  but  their  stronger  turn 
towards  the  surface  suggests  that  the  smaller  interaction  domain  results  in  stronger  turning  of 
the  outer  flow. 


Figure  ILI5.  Cross-stream  distributions  of  the  velocity  vector  at  x/c  =  0.16  and  0.21.  a=f  and 
Rcc  =  381,000.  fac  =  139Hz  and  850  Hz  (—). 

While  the  streamwise  position  of  the  suction  peak  appears  to  be  unaffected  by  the  actuation 
frequency,  its  magnitude  increases  with  the  frequency  and  decreases  with  the  jet  momentum 
coefficient.  The  flow  appears  to  be  most  sensitive  to  Cfi  stfact  =  139  Hz  (Figure  11. 14a)  where  a 
decrease  of  an  order  of  magnitude  in  the  momentum  coefficient  (to  1.7  10  )  results  in  flow 
separation  over  the  upper  side  of  the  airfoil  and  loss  of  the  actuation  effectiveness.  When  the 
actuation  frequency  is  increased  to  187  Hz  (Figure  II.  14b),  the  increase  in  the  suction  peak 
(Cn  =  -1 .95)  is  accompanied  by  a  decrease  in  the  streamwise  extent  of  the  recirculating  domain. 
When  Cn  is  decrease  by  an  order  of  magnitude  (to  1.7  10''’),  the  flow  separates  near  the  trailing 
edge  at  x/c  ~  0.8.  At  the  high  actuation  frequency  (850  Hz,  Figure  II.14d),  the  sensitivity  to  the 
level  of  Cfi  is  greatly  decreased  and  the  suction  peak  reaches  Cp  =  -2.4. 


Figure  11.16.  Pressure  drag  coefficient  (a)  and  lift  coefficient  (b)  as  a  function  of  the  actuation  frequency,  (f) 
Cp  =  1.2x10’/  (^)  5. 7x10’/  and  (*)  1. 7x10’/  Rec  =  381,000.  The  baseline  is  represented  by  the  dashed  line. 
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The  pressure  distributions  on  the  airfoil  at  the  various  actuation  frequencies  and  momentum 
coefficients  are  integrated  to  yield  the  lift  and  the  pressure  drag  coefficients  (Fibres  II.  16a  and 
b,  respectively,  where  the  baseline  levels  are  presented  by  the  dashed  lines).  It  is  interesting  to 
note  that  for  Cfi  =  1.2  10'^  the  pressure  drag  at  fact  -  139  Hz  actually  increases  by  a  factor  of  5.4 
relative  to  the  baseline  flow  as  a  result  of  the  streamwise  length  of  the  recirculating 
(interaction)  domain.  This  increase  becomes  significantly  larger  when  C,,  is  decreased. 
However,  as  the  actuation  frequency  is  increased,  for  a  given  momentum  coefficient  the 
pressure  drag  decreases  significantly  and  the  lift  increases.  These  data  also  show  that  for 
fact  =  850  Hz  the  actuation  leads  to  a  reduction  in  a  pressure  drag  and  some  lift  augmentation 
(except  for  C^=  1.7  10’'*,  where  the  lift  coefficient  is  unchanged)  and  the  trend  suggests  that 
further  increase  in  the  actuation  frequency  can  lead  to  additional  improvement  in  performance. 
As  noted  in  connection  with  Figure  11.14  above,  the  dependence  on  decreases  with 
increasing  actuation  frequency. 

The  effect  of  the  actuation  frequency  on  the  size  of  the  recirculating  (interaction)  domain  in  the 
vicinity  of  the  obstruction  is  investigated  using  particle  image  velocimetry  (PIV)  in  the  x-y 
plane  (z  =  0).  Each  measurement  field  is  comprised  of  four  partially  overlapping  frames 
measuring  10cm  x  10cm  (only  three  frames  are  used  for  /a«  =  850  Hz)  where  each  frame  is 
computed  from  an  ensemble  of  500  realization  pairs.  Figures  II.17a-d  show  magnified  views  of 
the  time-averaged  maps  of  cross  stream  vorticity  (that  also  include  regularly-spaced  cross 
stream  velocity  vectors  profiles)  within  the  domain  of  interaction  between  the  control  actuation 
and  the  cross  flow  for  fact~  139  Hz,  187  Hz,  362  Hz,  and  850  Hz,  respectively  (flr=3  , 
Rec  =  381,000  and  C^  =  1.2 10'^).  The  surface  of  the  airfoil  and  the  outline  of  the  obstruction 
are  included  in  each  frame  and  the  location  of  the  jet  which  is  placed  immediately  downstream 
of  the  obstruction  is  marked  with  an  arrow  [note  that  the  vertical  (y)  scale  is  magnified  1 .6: 1 
relative  to  the  horizontal  (x)  scale  to  show  more  detail].  In  the  presence  of  actuation,  the 
boundary  layer  upstream  of  the  obstruction  (having  a  characteristic  height  ft  =  0.01c  into  the 


Figure  11.17  Time-averaged  velocity  vector  and  spanwise  vorticity  maps  above  the  airfoil  at  a-  f, 
=  1.2x10-^  and  Rec  =  381,000 for  (a)f,c,  =  139  Hz,  (b)  187,  (c)  362  Hz,  and  (d)  850  Hz. 
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cross  flow)  separates  locally  and  a  closed  recirculation  domain  is  formed  as  a  result  of  the  jet 
actuation.  The  images  show  the  clockwise  (CW)  spanwise  vorticity  concentration  within  the 
boundary  layer  and  the  ensuing  separating  shear  layer,  as  well  as  the  counter  clockwise  (CCW) 
vorticity  near  the  surface  within  the  recirculation  domain  downstream  of  the  obstruction.  As 
noted  above,  the  streamwise  and  cross  stream  extent  of  this  domain  are  influenced  by  the 
frequency  and  strength  of  the  jet.  Actuation  at  139  Hz  (Figure  II.  17a)  results  in  a  recirculating 
domain  that  extends  almost  to  the  downstream  edge  of  the  frame  or  26  h  downstream  of  the 
obstruction  and  2.2  h  into  the  cross-stream  direction.  Increasing  the  actuation  frequency  to  187 
Hz  (Figure  II.  17b)  results  in  a  somewhat  smaller  interaction  domain  having  typical  streamwise 
len^h  and  cross-stream  height  of  20  /j  and  1.7 /i,  respectively.  When  the  actuation  frequency  is 
increased  to  362  Hz  (Figure  II.  17c)  and  850  Hz  (Figure  II.17d)  it  decreases  the  respective 
characteristic  streamwise  and  cross-stream  lengths  of  the  recirculating  domain  to  147i  and  I  h 
and  1.3 and  \  h.  As  shown  by  the  vector  maps,  the  smaller  interaction  domain  (Figure  II.  17c) 
leads  to  a  stronger  turning  of  the  outer  flow  and  therefore  results  in  a  stronger  and  sharper  low- 
pressure  (suction)  peak  on  the  surface  of  the  airfoil. 

The  nominal  streamwise  length  of  the  recirculating  domain  Z,  is  estimated  based  on  the 
reattachment  location  of  the  recirculating  domain  downstream  of  the  obstruction  and  is 
normalized  by  the  chord  of  the  airfoil.  As  shown  in  Figure  11.18,  LJc  decay  exponentially  with 
F^.  It  is  noteworthy  that  when  the  jet  momentum  coefficient  is  reduced  by  (roughly)  an  order 
of  magnitude  (to  1.710'^),  its  effectiveness  at  low  actuation  frequencies  diminishes 
dramatically.  This  is  demonstrated  in  a  sequence  of  images  in  Figures  II.19a-c.  While  for 
/art  =  139  Hz  the  flow  above  the  airfoil  is  separated  (not  shown),  at  187  Hz  the  recirculation 
domain  extends  through  0.8c  (based  on  the  pressure  measurements  in  Figure  II.  14b)  which  is 
nominally  2  ?iact  wavelengths  of  this  actuation  frequency  (for  which  =  1.7).  However,  when 


Figure  11.18.  Normalized  streamwise  length  of  the  recirculating  domain  as  a  Junction  of  the  non- 
dimensional Jrequency.  a  =  /  and  Rcc  =  381,000. 

C^=1.210''’  (Figure  II.  17b),  the  length  of  the  circulation  domain  for  the  same  actuation 
frequency  is  only  0.58  Xact  (cf.  Figure  11.18)  indicating  that  the  scaling  of  the  interaction  domain 
is  not  merely  a  function  of  the  actuation  frequency  and  depends  also  on  the  jet  momentum 
coefficient.  Similarly,  for  /rt=362  Hz,  the  streamvvise  length  of  the  circulation  domain 
increases  from  0.82  X.act  to  1.8  and  for  fact  =  850  Hz  it  increases  from  1.1  X^ct  to  2Xact  when 
is  decreased  (i.e.,  compare  Figures  II.  17c  and  19b  and  17d  and  19c). 
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The  structure  of  the  interaction  domain  for  /acr  =  850  Hz  is  investigated  using  phase  averaged 
PIV  measurements  (phase-locked  to  the  actuation  wavefom).  The  images  in  Figures  II.20a-f 
are  acquired  at  six  equal  time  increments  during  the  actuation  cycle  where  Figures  II.20a  and  d 
correspond  to  the  beginning  of  the  blowing  and  suction  strokes,  respectively.  At  the  beginning 
of  the  blowing  stroke  (Figure  IL20a),  the  recirculation  domain  includes  a  CW  vortex  centered 
approximately  at  x'=  10  h  downstream  from  the  obstruction.  A  smaller  CW  vortex  is  also 
visible  farther  downstream  (x'  =  28 Ti).  As  the  blowing  proceeds  (Figures  II.20b  and  c),  the 


Figure  11.19.  Time-averaged  velocity  vectors  and  spartwise  vorticity  maps  above  the  airfoil  at 
a  =  3°,  Q  =  1. 7x10-'  and  Re,  =  381,000 for  (a)fac,  =  187  Hz.  (b)  362,  and  (c)  850  Hz. 

upstream  (larger)  CW  vortex  seems  to  intensify  as  is  evidenced  by  the  relatively  strong 
downwash  that  is  induced  near  its  downstream  edge  and  it  is  advected  farther  downstream  (the 
maximum  streamwise  length  of  the  primary  recirculating  region  is  ~7 7?).  Note  the  formation  a 
counter  rotating  vortex  pair  during  the  blowing  cycle  just  downstream  of  the  obstruction. 
Figure  II.20b).  At  the  same  time,  it  appears  that  the  strength  of  the  secondary  (downstream) 
vortex  is  somewhat  reduced.  The  effect  of  the  suction  stroke  is  apparently  to  ingest  the  weak 
CCW  vortex  (immediately  downstream  of  the  obstruction)  and  to  induce  the  formation  of  a 
second  CW  vortex  within  the  primary  recirculating  domain  (Figure  II.20e),  which  intensifies  as 
the  downstream  (earlier)  CW  vortex  becomes  weaker.  By  the  time  the  cycle  repeats  itself,  the 
latter  vortex  is  advected  outside  of  the  downstream  edge  of  the  primary  recirculating  domain. 
This  sequence  of  images  clearly  indicates  that  for  this  level  of  the  streamwise  length  of  the 


25 


interaction  domain  is  of  the  order  of  one  to  two  wavelengths  of  the  actuation  frequency.  More 
importantly,  similar  to  the  earlier  observations  of  Honohan  et  al.'^  farther  downstream,  the 
actuation  frequency  has  negligible  temporal  effect  on  the  vorticity  distribution  within  the  wall 
boundary  layer. 

In  order  to  further  explore  the  effect  of  the  actuation  frequency  on  the  formation  of  the 
interaction  domain,  phase-averaged  images  similar  to  the  data  presented  in  Figure  11.20  are 
shown  in  Figures  II.21a-d  for  fact  -  139  Hz,  187  Hz,  362  Hz,  and  850  Hz,  respectively.  These 
phase-averaged  images  are  acquired  at  the  beginning  of  the  blowing  stroke  for  each  actuation 
frequency.  For  139  Hz  (Figure  11.21  a)  the  recirculation  domain  includes  a  large  CW 
vortex  centered  approximately  at  jc'  ~  33  downstream  from  the  obstruction  (a  second  weaker 
vortex  CW  vortex  is  present  downstream  of  the  measurement  domain).  PA  fact  =  187  Hz 
(Figure  II.21b),  the  time-averaged  recirculation  domain  in  Figure  II.  17b  includes  two  CW 
vortices  centered  approximately  at  x'=  107?  and  247?  downstream  from  the  obstruction,  where 
the  downstream  vortex  appears  to  be  considerably  more  diffused  into  the  free  stream 
(ostensibly  owing  to  loss  of  phase  coherence)  and  is  closer  to  the  downstream  edge  of  the 
domain  outlined  by  the  time-averaged  data  in  Figure  II.  17b.  Further  increase  of  the  actuation 
frequency  to  362  Hz  and  850  HZ  (Figures  11.21  c  and  d,  respectively)  results  in  smaller  vortices 
that  are  located  at  x'  =  \0  h  and  20  h,  and  x'  =  5  h,  and  \2  h,  respectively.  A  comparison  of  the 


Figure  11.20.  Phase-averaged  velocity  vectors  and  spanwise  maps  above  the  airfoil  at  a  =  i®,  =  1.2x1 

u,  =  850  Hz  and  Re,  =  381,000.  (a)  <l>  =  0P.(b)  6(f,  (c)  120‘>,  (d)  18(f.  (e)  24(f,  and  (f)  30(f. 
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images  in  Figure  11.21  with  the  corresponding  time-averaged  images  in  Figure  11.17  indicates 
that  the  bulk  of  the  interaction  domain  scales  approximately  with  one  to  two  vortical  stnictures 
of  the  actuation  wavelength  where  the  earlier  (downstream)  vortex  becomes  considerably 
weaker  and  loses  its  identity  as  it  is  advected  downstream.  It  is  also  noteworthy  that  at  low 
actuation  frequencies  the  cross-stream  extent  of  the  interaction  domain  is  larger  than  the 
obstruction  height,  which  might  explain  the  blockage  (velocity  reduction)  that  is  evident  from 
the  pressure  data  in  Figure  11. 14. 
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Figure  11.21  Phase-averaged  velocity  vectors  and  spanwise  vorticity  maps  above  the  airfoil  at  a- 3^, 
=  IJxIOr^  and  Re,  =  381, 000 for  (a)fa„  =  139  Hz.  (b)  187  Hz.  (c)  362  Hz  and  (d)  850  Hz. 

11.6.  CONCLUSIONS 

The  objective  of  the  control  approach  described  here  is  to  reduce  the  pressure  drag  at  low 
angles  of  attack  with  minimal  or  no  penalty  in  lift.  The  pressure  distribution  on  the  surface  of 
the  airfoil  is  modified  using  virtual  aero-shaping  that  is  achieved  by  the  formation  of  a 
stationary  recirculating  flow  domain  next  to  the  surface  as  a  result  of  the  interaction  be^een 
the  cross  flow  and  a  high-frequency  synthetic  jet  actuator  placed  downstream  from  a  miniature 
surface-mounted  passive  obstruction.  The  formation  of  this  interaction  domain  alters  the  flow 
above  the  airfoil  by  inducing  a  displacement  of  the  local  streamlines  that  is  sufficient  to  modify 
the  surface  pressure  distribution  even  at  small  angles  of  attack. 

The  present  (and  earlier)  measurements  show  that  while  relatively  low-level  actuation  [i.e., 
of  0(10'^)]  can  effectively  induce  fully  attached  flow  at  post  stall  angles  of  attack  of  the 
baseline  airfoil,  the  induced  aerodynamic  changes  at  low  angles  of  attack  ^e,  in  effect, 
marginal.  As  shown  in  earlier  work,  control  effectiveness  can  be  improved  by  increasing  the 
actuation  levels  (i.e.,  higher  Cn)  as  may  be  achieved  by  arrays  of  jet  actuators.  However,  the 
present  approach  is  to  place  a  low-C^  synthetic  jet  actuator  downstream  from  a  miniature 
surface-mounted  passive  obstruction  that  typically  extends  0.01c  above  the  surface  of  the 
airfoil.  Measurements  of  the  pressure  distribution  about  the  surface  of  the  airfoil  show  that 
compared  to  the  baseline  configuration  actuation  within  the  domain  0.19  <  x/c  <  0.28  results  in 
increased  suction  in  the  vicinity  of  the  actuator  and  a  reduction  in  the  magnitude  of  the  pressure 
within  the  recovery  domain  towards  the  trailing  edge.  These  changes  result  in  reductions  in  the 
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pressure  drag  that  are  comparable  to  the  magnitude  of  the  pressure  drag  of  the  baseline 
configuration,  and  concomitantly  in  relatively  small  reductions  (between  nil  and  8%)  in  lift. 

In  the  present  work,  the  height  of  the  obstruction  (above  the  surface  of  the  airfoil)  is 
h  =  9.310'^c  (i.e.,  less  than  1%  of  the  chord)  and  it  is  placed  at  Xc/c  =  0.22  where  the  synthetic 
jet  actuator  is  located  immediately  downstream  of  the  obstruction.  The  experiments  are 
conducted  at  Rec  =  381,000  for  which  the  baseline  flow  is  fully  attached  to  the  surface  of  the 
airfoil  for  a<  8°. 

The  interaction  domain  between  the  actuator  and  the  cross  flow  is  investigated  using  high- 
magnification  particle  image  velocimetry  in  the  x-y  plane  (z  =  0).  The  characteristic  streamwise 
scale  of  the  induced  recirculating  flow  domain  is  varied  by  changing  the  actuation  frequency 
(F^  is  varied  between  1.25  and  7.8)  and/or  the  jet  strength  (Q/  varies  between  1.7 10  and 
1.2 10'^).  The  present  work  shows  that  the  streamwise  length  of  the  interaction  domain  and  its 
cross  stream  penetration  increase  with  decreasing  actuation  frequency  and  jet  momentum 
coefficient.  Phase-averaged  PIV  images  show  that  for  Q/  =  1.2 10' ,  at  given  actuation 
frequency  the  streamwise  length  of  the  interaction  domain  is  typically  comparable  to  one  to  two 
wavelengths  of  the  actuation  frequency.  The  interaction  domain  contains  one  strong  clockwise 
vortex  and  a  second  weaker  vortex  (from  the  earlier  actuation  cycle)  is  advected  downstream 
and  rapidly  loses  its  identity.  However,  a  reduction  in  C^i  by  an  order  of  magnitude  (to  1.7 10  ) 
significantly  lengthens  the  interaction  domain  indicating  that  it  does  not  merely  scale  with  the 
actuation  frequency.  For  example,  at  F*  =  7.8  (850  Hz)  the  length  of  the  interaction  domain 
increases  from  XA^ct  to  iXact  and  at  lower  actuation  frequency  for  which  F^  is  1.25  the 
decrease  in  Cfi  results  in  a  mismatch  between  the  interaction  domain  and  flow  separation. 

As  a  result  of  the  changes  in  the  curvature  of  the  interaction  domain,  the  induced  low-pressure 
(suction)  peak  on  the  airfoil  diminishes  and  widens  with  the  actuation  frequency  and  jet 
momentum  coefficient  leading  to  a  substantial  increase  in  pressure  drag  and  a  reduction  in  lift. 
It  is  remarkable  that  a  significant  decrease  in  the  baseline  pressure  drag  (compared  to  its 
magnitude)  with  minimal  penalty  in  lift  is  only  achieved  at  higher  actuation  frequencies 
^c/  =  850  Hz).  The  present  results  also  show  that  at  high  actuation  frequencies  the 
aerodynamic  forces  (and  the  shape  of  the  interaction  domain)  are  substantially  less  sensitive  to 
changes  in  the  jet  momentum  coefficient. 
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III.  AERODYNAMIC  MODIFICATION  AT  LOW  ANGLES  OF  ATTACK  USING 
TRAILING  EDGE  TRAPPED  VORTICES 

III.1.  OVERVIEW 

The  aerodynamic  performance  of  a  conventional  non-symmetric  airfoil  is  manipulated  at  low 
angles  of  attack  (when  the  baseline  flow  is  fully  attached)  by  leveraging  the  presence  of  a  small 
flap  integrated  with  a  synthetic  jet  actuator.  The  flap  is  similar  in  some  respects  to  a 
conventional  Gurney  flap  having  a  characteristic  scale  of  0.016c  and  is  mounted  on  the  pressure 
side  of  the  airfoil  just  upstream  of  its  trailing  edge.  The  objective  is  to  leverage  the  presence  of 
the  trapped  vortices  that  are  induced  by  the  flap  near  the  trailing  edge  to  further  manipulate  the 
Kutta  condition  of  the  airfoil  and  thereby  achieve  global  flow  modifications  and  potentially 
improve  it's  overall  aerodynamic  performance.  When  the  control  jet  is  activated,  the 
aerodynamic  performance  of  the  airfoil  can  be  continuously  varied  from  a  configuration  that  is 
essentially  an  airfoil  with  a  conventional  Gurney  flap  through  an  airfoil  that  has  similar 
performance  as  that  of  the  smooth  (baseline)  airfoil  to  an  airfoil  that  has  a  higher  lift  to  pressure 
drag  ratio. 

III.2  INTRODUCTION 

Utilization  of  flow  control  techniques  in  aerodynamic  applications  has  traditionally  concentrated 
on  separated  flows  over  stalled  wing  sections  or  flaps  with  the  objective  of  improving  the 
aerodynamic  performance  by  inducing  some  degree  of  flow  attachment.  Control  approaches 
have  relied  either  on  exploiting  the  instability  of  the  separating  shear  layer  at  relatively  low 
[F^~0(1)]  actuation  frequencies  (e.g.,  Hsiao  et  al.*,  Seifert  et  al.^,  Chang  et  al.^),  or  the 
suppression  of  separation  by  an  "apparent"  modification  of  the  flow  boundary  at  high 
[F"^  ~  0(10)]  actuation  frequencies  (e.g.,  Smith  et  ah'*  Amitay  et  al.^’^  and  Erk^). 

The  recent  works  of  Chatlynne  et  al.®  and  Amitay  et  al.^  has  demonstrated  that  the  modification 
of  aerosurface  by  means  of  high-frequency  fluidic  actuation  can  be  utilized  for  aerodynamic  flow 
control  even  in  the  absence  of  stall.  This  is  particularly  important  for  cruise  conditions  at  low 
angles  of  attack  where  the  baseline  flow  is  mostly  attached.  These  authors  achieved  drag 
reduction  with  minimal  lift  penalty  on  a  Clark-Y  airfoil  at  low  angles  of  attack.  The  flow  above 
the  airfoil  (i.e.,  on  the  suction  side)  was  altered  by  the  formation  of  a  stationary  recirculating  flow 
domain  near  the  surface  using  a  high-frequency  synthetic  jet  actuation  downstream  from  a 
miniature  (nominally  0.01c)  surface-mounted  passive  obstruction.  The  formation  of  a  trapped 
vortex  within  the  Interaction  domain  between  the  synthetic  jet  and  the  obstruction  results  in 
reduction  in  the  pressure  drag  that  is  comparable  to  the  magnitude  of  the  pressure  drag  of  the 
baseline  conflguration  and  is  accompanied  by  a  relatively  small  reduction  (between  nil  and  8%) 
in  lift. 

This  work  is  a  continuation  of  the  earlier  work  on  aerodynamic  flow  control  at  low  angles  of 
attack  (§11).  While  the  work  of  Chatlynne  et  al.*  and  Amitay  et  al.^  focused  on  flow  modiflcation 
on  the  pressure  side  of  the  airfoil  near  the  leading  edge,  the  current  work  exploits  the  presence  of 
a  miniature  Jet  Integrated  Flap  (JIF)  similar  to  a  Gurney  flap  that  is  placed  on  the  pressure  side  of 
the  airfoil  just  upstream  of  the  trailing  edge. 
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III.3  THE  TRAILING  EDGE  CONTROL  CONCEPT 

The  Gumey  flap  is  a  flat  plate  that  is  typically  mounted  along  the  trailing  edge  of  a  conventional 
airfoil  and  normal  to  its  lower  (pressure)  surface  such  that  its  height  is  nominally  1  to  2  percent 
of  the  chord.  As  noted  by  Liebeck'®  (and  confirmed  by  wind  tunnel  testing  of  a  flap  of  height 
0.0125c  on  the  symmetric  Newman  airfoil),  when  the  flap  height  is  less  than  0.02c,  it  can  lead  to 
simultaneous  increase  in  lift  and  reduction  in  drag.  Although  Liebeck'®  noted  that  the 
mechanism  for  drag  reduction  was  (and  apparently  still  is)  not  well  understood,  he  nevertheless 
suggested  a  possible  mechanism  as  outlined  schematically  in  Figures  III. la  and  b.  Liebeck’" 
conjectured  that  the  drag  reduction  results  from  the  displacement  of  small  adjacent  separation 
bubbles  (essentially  a  counter-rotating  vortex  pair)  that  form  along  the  upper  and  lower  surfaces 
of  the  airfoil  just  upstream  of  its  trailing  edge  (Figure  III.  la).  These  bubbles  extend  into  the  near 
wake  and  a  stagnation  point  is  formed  along  the  dividing  streamline  that  extends  to  the  trailing 
edge.  The  downwash  associated  with  the  increased  lift  in  the  presence  of  the  Gumey  flap 
apparently  results  in  displacement  and  stretching  of  the  separation  bubbles  (Figure  III. lb)  and 
therefore  in  a  smaller  velocity  deficit  within  the  airfoil's  wake.  The  formation  of  the  counter 
rotating  vortex  pair  downstream  of  a  Gumey  flap  was  later  demonstrated  computationally  on  a 
cambered  NACA  4412  airfoil  profile  by  Jang"  and  confirmed  by  LDV  measurements  of  Jeffrey 
et  al.*^,  who  also  show  how  a  Gumey  flap  increases  circulation  over  the  airfoil  by  functioning  as 
the  equivalent  of  a  stationary  vortex  located  at  the  trailing  edge,  and  give  further  insight  into  the 
vortex  shedding  mechanism.  Giguere  et  al’^  conjectured  that  drag  reduction  can  be  attained  if 
the  height  of  the  Gumey  flap  is  smaller  than  the  local  boundary  layer  thickness.  There  is  also 
some  evidence  in  the  literature  (e.g.,  Myose  et  al.'"*)  that  suggests  that  drag  reduction  by  Gumey 
flaps  is  more  readily  achieved  on  symmetric  airfoils. 


Counter-Rotating 
Vortex  Pair 


Stagnation  Point 


Figure  HI.l.  Hypothesized  trailing-edge  flow  conditions  of  an  airfoil  with  a  Gurney  flap.  (From  Liebeck!'’). 


The  objective  of  the  present  work  is  to  manipulate  the  aerodynamic  performance  of  a 
conventional  airfoil  by  leveraging  the  presence  of  a  Gumey  flap  near  its  trailing  edge. 
Specifically,  the  idea  is  to  manipulate  the  trapped  vortices  that  are  induced  by  the  flap  to  alter  the 
Kutta  condition  and  achieve  global  flow  modifications.  As  noted  in  Section  III.2,  the  concept  of 
using  trapped  vortices  to  alter  the  aerodynamic  performance  of  aero-surfaces  is  not  new  and  was 
already  addressed  in  the  1940s  and  50s.  One  example  of  this  approach  is  the  use  of  a  stationary. 
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trapped  vortex  to  alter  the  apparent  local  surface  curvature  and  therefore  the  direction  of  the 
external  flow.  Perkins  and  Hazen'^  used  a  suction  slot  near  the  trailing  edge  of  an  airfoil  and  an 
upstream  cavity  to  trap  a  vortex  near  the  trailing  edge  and  reported  substantial  lift  {Ci  =  1 .2)  at 
zero  angle  of  attack  at  relatively  low  (3%)  suction  coefficient.  Mandl’^  used  a  split  flap  located 
on  the  pressure  side  of  an  airfoil  combined  with  suction  (at  the  flap)  to  trap  a  vortex  and  to 
enforce  tangential  flow  near  the  trailing  edges  of  the  airfoil  and  the  flap  and  alter  the  Kutta 
condition  of  the  airfoil. 

III.4  EXPERIMENTAL  APPARATUS  AND  PROCEDURE 

The  present  experiments  use  a  modular  NACA  4415  airfoil  model.  The  2-D  airfoil  (shown 
schematically  in  Figure  III.2)  has  a  chord  measuring  461  mm  and  maximum  thickness  to  chord 
ratio  t/c  =  0.15  and  spans  the  entire  wind  tunnel  test  section  (which  measures  approximately  1  m 
on  the  side).  The  center  section  of  airfoil  model  is  equipped  with  two  configurations  of  modular 
spanwise  synthetic  jet  actuators  each  having  a  rectangular  orifice  measuring  140  x  0.5  mm.  One 
configuration  of  actuators  can  be  installed  within  the  streamwise  domain  0.09  <x/c  <0.315 
where  the  Jets  are  flush  mounted  and  issue  normal  to  the  surface  of  the  airfoil.  In  a  second 
configuration,  a  separate  jet  actuator  is  integrated  with  a  miniature  flap  having  a  characteristic 
height  of  0.016c.  The  Jet  Integrated  Flap  (JIF)  is  movable  along  the  surface  of  the  airfoil  on  the 
pressure  side  (Figure  III.2).  In  both  configurations,  the  actuator  jets  are  driven  by  several 
piezoelectric  membranes  that  are  built  into  a  central  cavity  and  are  operated  at  or  near  resonance. 
For  both  actuator  configurations,  streamwise  fences  are  placed  near  the  (spanwise)  edges  of  the 
actuator  orifices  to  maintain  a  nominally  two-dimensional  flow  in  the  controlled  section  of  the 
airfoil. 


y 


c  =  457  mm 


Control  Jet 

Figure  JTL2.  Airfoil  model 

The  airfoil  is  also  instrumented  with  a  circumferential  array  of  71  pressure  ports  (located  at  mid¬ 
span)  that  are  each  connected  to  an  external  high-speed  pressure  measurement  system. 

The  velocity  and  vorticity  fields  in  the  cross-stream  (x-^-)  plane,  z  =  0,  above  the  airfoil  (i.e.,  on 
the  suction  side)  and  in  its  wake  are  measured  using  particle  image  velocimetry  (PIV).  The  flow 
is  seeded  with  micron  smoke  particles  and  it  is  illuminated  using  a  double-pulse  ND-YAG  laser. 
Image  pairs  are  captured  using  a  1008  x  1016  pixel  CCD  camera  with  a  magnification  of  55.2 
pm/pixel.  In  each  image,  the  velocity  vectors  are  computed  using  a  standard  cross-correlation 
technique.  A  given  PIV  map  is  averaged  over  500  realizations  (image  pairs). 

Iil.5  AERODYNAMIC  PERFORMANCE  WITH  ACTIVE  GURNEY  FLAP 

As  discussed  in  Section  III.3,  the  addition  of  a  Gurney  flap  that  scales  on  the  order  of  a  few 
percent  of  the  chord  length  of  a  conventional  airfoil  can  lead  to  an  increase  in  lift  that  can  be 
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accompanied  by  an  increase  or  decrease  in  pressure  drag.  In  the  present  work,  the  JIF  is  similar 
in  some  respects  to  a  conventional  Gurney  flap  and  is  mounted  on  a  conventional  (smooth) 
airfoil.  The  objective  is  to  leverage  the  presence  of  the  obstruction  (and  the  altered  baseline  flow 
of  the  smooth  airfoil)  to  further  manipulate  the  Kutta  condition  of  the  airfoil  to  achieve  global 
flow  modifications  and  potentially  improve  the  airfoil's  overall  aerodynamic  performance.  An 
important  element  of  this  approach  is  that  it  is  implemented  at  low  angles  of  attack  and  requires 
relatively  low  actuation  input  by  leveraging  the  presence  of  the  flap  as  a  baseline. 

The  aerodynamic  characteristics  of  the  present  JIF  was  compared  with  a  "conventional"  0.014c 
surface-normal  Gurney  flap  when  both  are  mounted  flush  with  the  trailing  edge  of  the  baseline 
(smooth)  airfoil.  The  bulk  of  the  experiments  reported  here  are  conducted  at  Rcc  =  460,000 
(based  on  the  free  stream  velocity  and  the  chord  of  the  airfoil).  Figures  III.3a-c  show  the 
variations  with  angle  of  attack  (-4°  <  flr<  8°)  of  the  lift  and  (pressure)  drag  coefficients  {Cl  and 
Cop)  and  of  the  lift  to  pressure  drag  ratio  {UDp),  respectively,  while  Figure  111.3d  is  the 
corresponding  pressure  drag  polar.  Data  for  the  baseline  (smooth)  airfoil,  the  airfoil  combined 
with  the  conventional  Gurney  flap,  and  the  JIF  are  shown.  The  addition  of  a  conventional 
Gurney  flap  leads  to  an  increase  in  the  lift  coefficient  (nominally  0.12)  that  is  almost  invariant 
over  the  entire  range  of  angles  of  attack  tested  here  (Figure  III.3a).  As  shown  in  Figure3b,  the 
addition  of  a  Gurney  flap  to  the  present  (cambered)  airfoil  results  in  an  increase  in  pressure  drag 


Figure  III.3.  The  variations  with  angle  of  attack  of  the  Cl  (a),  Cop  (b),  and  L/Dp  (c)  and  the  corresponding 
pressure  drag  polar  (d).  (O)  Smooth,  (A)  Gurney  flap  and  (O)  JIF. 
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from  112%  at  cr=0°  to  46%  at  a=8°.  Similar  trends  were  also  observed  in  the  studies  on 
cambered  airfoils  by  Jeffrey  et  al.*^  and  Giguere  et  al.'^  as  well  as  in  the  CFD  results  of  J  Jang  . 
In  these  cases,  Cl  and  Co  increase  at  a  given  angle  of  attack  with  addition  of  the  Gurney  flap. 

Despite  the  increase  in  lift,  the  addition  of  the  Gurney  flap  at  the  trailing  edge  results  in  a 
decrease  in  the  lift  to  drag  ratio  (L/Dp)  for  a>  -2°.  For  smaller  angles  of  attack,  UDp  with  the 
Gurney  flap  is  higher  than  for  the  smooth  airfoil.  The  data  in  Figure  111.3  show  that  the 
aerodynamic  characteristics  of  the  baseline  airfoil  combined  with  the  present  jet  integrated  flap 
are  reasonably  similar  to  what  is  measured  with  the  conventional  Gurney  flap.  While  the  JIF 
yields  a  slightly  higher  lift  coefficient  for  a>  2°  (Figure  III.3a),  there  is  also  a  small  increase  in 
the  drag  coefficient  (Figure  II1.3b)  such  that  L/Dp  is  almost  the  same  within  the  range  2°  <  a  <  8®. 

The  effect  of  jet  actuation  on  the  pressure  distribution  around  the  airfoil  when  the  JIF  is  placed  at 
xq/c  =  0.94  is  shown  in  Figures  111.4a  and  b  for  a=  0®  and  8®,  respectively,  (the  location  of  the 
downstream  edge  of  the  JIF  is  marked  with  a  dashed  line).  The  pressure  distribution  in  the 
absence  of  actuation  (open  symbols)  shows  that  for  x/c  <  0.8  and  similar  to  a  conventional 
Gurney  flap  as  indicated  by  Giguere  et  al.'^,  the  unactuated  JIF  leads  to  respective  decrease  and 
increase  in  the  pressure  coefficient  (which  corresponds  to  an  increase  in  the  lift)  on  the  suction 
and  pressure  surfaces  of  the  airfoil.  It  is  remarkable  that  these  changes  are  measurable  as  far 
upstream  as  x/c  =  0.1.  Compared  to  the  smooth  airfoil,  the  increase  in  lift  that  results  from  the 
decrease  in  pressure  downstream  from  the  suction  peak  on  the  top  surface  is  ostensibly 
accompanied  by  increased  (downward)  turning  of  the  flow  above  the  airfoil  and  a  concomitant 


Figure  III.4.  The  pressure  distribution  around  the  airfoil  for  a  =  (T  (a)  and  8°  (b).  The  downstream  edge  of  the  JIF 
is  marked  with  a  dashed  line  (xq/c  =  0. 94).  ( - )  Smooth,  (tT)  unactuated  JIF  and  (  )  actuated  JIF. 
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reduction  in  the  pressure  coefficient  at  the  trailing  edge.  However,  as  noted  in  connection  with 
Figure  IIL3b,  the  increase  in  lift  is  also  accompanied  by  an  increase  in  pressure  drag.  The 
presence  of  the  JIF  on  the  pressure  surface  of  the  airfoil  results  in  slight  blockage  that  peaks  near 
x/c  =  0.80  followed  by  an  almost  linear  decrease  of  the  pressure  towards  the  trailing  edge. 

The  actuation  of  the  jet  (C;/=  8x10'^)  results  in  a  localized  low  pressure  domain  between  the 
downstream  edge  of  the  actuator  and  the  trailing  edge  of  the  airfoil  as  is  evidenced  by  the 
presence  of  a  suction  peak  at  x/c  =  0.93  where  the  pressure  coefficient  is  lower  than  in  the 
absence  of  actuation  while  the  pressure  at  the  trailing  edge  of  the  airfoil  is  higher  when  the  jet  is 
active.  In  fact,  it  is  remarkable  that  when  the  jet  is  active,  the  pressure  distribution  along  the 
entire  upper  (suction)  surface  of  the  airfoil  is  lower  than  in  the  absence  of  actuation  and  is 
virtually  identical  to  the  pressure  distribution  of  the  smooth  airfoil.  At  the  same  time,  while  the 
pressure  distribution  along  the  pressure  side  of  the  airfoil  is  generally  lower  than  when  the  jet  is 
inactive,  it  is  still  higher  than  for  the  smooth  airfoil.  Up  to  x/c  =  0.80  (where  there  is  a  local 
maximum)  the  pressure  distributions  in  the  presence  and  absence  of  the  actuation  are  rather 
similar.  However,  the  low-pressure  domain  that  is  present  near  the  airfoil's  trailing  edge  appears 
to  migrate  upstream  towards  the  JIF  when  the  jet  is  activated.  As  will  be  shown  in  Figure  III.9 
below,  the  actuation  results  in  upward  turning  of  both  flows  along  the  lower  (pressure)  and  upper 
(suction)  surfaces  towards  the  trailing  edge,  which  accounts  for  the  increase  in  the  trailing  edge 
pressure.  The  upward  tilt  is  clearly  commensurate  with  the  reduction  in  lift  that  is  caused  by  the 
appearance  of  the  suction  peak  on  the  lower  surface  {x/c  =  0.93).  The  upward  turning  of  the  flow 
near  the  trailing  edge  is  similar  to  the  turning  (vectoring)  of  a  primary  jet  towards  a  low-pressure 
interaction  domain  that  is  induced  by  an  adjacent  synthetic  jet  actuator  which  was  investigated  in 
detail  by  Smith  and  Glezer'^. 


Figure  III.5.  (a)  The  variation  ofUDp  with  angles  of  attack,  and  (b)  the  pressure  drag  polar.  (O)  Smooth,  ( O) 
unactuated  JIF  and  ( 0)  actuated  JIF. 


The  effect  of  the  actuation  on  the  aerodynamic  characteristics  of  the  airfoil  for  angles  of  attack 
within  the  range  -2°  <  a<  8”  is  shown  in  Figure  III.5.  As  shown  in  Figure  III.3c,  although  the 
addition  of  the  JIF  at  xo/c  =  0.94  results  (in  the  absence  of  actuation)  in  an  increase  in  lift, 
however,  an  increase  in  pressure  drag  leads  to  an  overall  decrease  in  UDp  compared  to  the 
baseline  (smooth)  airfoil.  When  actuation  is  applied  (Q  =  8x10'^),  there  is  a  substantial  increase 
in  L/Dp  (Figure  Ill.Sa)  compared  to  both  the  baseline  and  the  unactuated  JIF  airfoils  over  the 
entire  range  of  angles  of  attack.  At  this  streamwise  position  of  the  JIF,  the  maximum  lift  to 
pressure  drag  ratio  is  obtained  at  Of=  0“  {L/Dp  =  36.3),  which  is  an  increase  of  45%  over  the 
unactuated  JIF  airfoil,  but  more  significantly  an  increase  of  15.3%  over  the  smooth  airfoil.  The 
corresponding  pressure  drag  polar  (Figure  III.5b)  shows  that  for  a  given  lift  coefficient,  there  is  a 
decrease  in  the  pressure  drag  compared  to  both  the  unactuated  JIF  and  smooth  airfoils.  As  Cl 
increases,  the  decrement  in  pressure  drag  increases  relative  to  the  smooth  airfoil  and  decreases 
relative  to  the  unactuated  JIF  airfoil. 


Figure  III.6.  Variation  with  Q  of  the  normalized  increments  in  the  (a)  lift,  (b)  pressure  drag,  and  (c)  L/D. 


36 


The  effectiveness  of  the  actuation  can  be  continuously  varied  with  the  magnitude  of  the  control 
jet  momentum  coefficient  (e.g.,  Amitay  et  al.^).  Figures  III.6a-c  show  the  variation  with  of 
the  normalized  increments  (relative  to  the  smooth  airfoil)  in  the  lift,  pressure  drag,  and  lift  to 
drag  ratio,  respectively  at  a=  0°.  For  the  unactuated  JIF  case  (i.e.,  =  0),  the  respective 

increases  in  lift  and  drag  are  35%  and  68%,  and  the  corresponding  decrease  in  L/Dp  is  20%.  As 
Cji  is  increased,  there  is  a  monotonic  decrease  in  both  the  normalized  lift  and  pressure  drag  and 
an  increase  in  L/Dp.  Apparently,  the  rate  of  decrease  with  Cp  of  the  normalized  pressure  drag  is 
larger  than  that  of  the  normalized  lift,  and  the  lift  to  pressure  drag  ratio  of  the  smooth  airfoil  is 
attained  at  Q  =  4.9x10'^  when  the  lift  and  drag  coefficients  are  each  8%  larger  than  for  the 
smooth  airfoil.  The  increase  in  pressure  drag  for  the  actuated  JIF  case  vanishes  at  Q/*  6x10' 
while  Cl  is  still  4.5%  higher  than  for  the  smooth  airfoil.  Compared  to  the  smooth  airfoil,  at  the 
highest  actuation  level  that  was  achieved  with  the  present  actuators  (C^=  8x10'^),  L/Dp  increases 
by  1 9%,  the  pressure  drag  decreases  by  1 6%  while  the  lift  is  nominally  unchanged. 

These  data  show  that  the  aerodynamic  performance  of  a  conventional  airfoil  can  be  continuously 
varied  by  adding  a  small  (Gurney  flap  like)  obstruction  integrated  with  a  fluidic  actuator.  In 
fact,  the  aerodynamic  performance  of  the  airfoil  can  be  varied  from  a  configuration  that  is 
essentially  an  airfoil  with  a  conventional  Gurney  flap  through  an  airfoil  that  has  similar 
performance  as  that  of  the  smooth  (baseline)  airfoil  to  an  airfoil  that  has  higher  lift  to  pressure 
drag  ratio. 

The  effects  of  the  actuation  on  the  airfoil's  moment  coefficient  and  on  the  center  of  pressure 
(about  c/4)  are  shown  in  Figures  Ilf.Ta  and  b,  respectively  for  -2‘’<£ir<8°.  The  center  of 


r^ure  in.  7.  The  variation  with  a  of  the  (a)  moment  coefficient  and  (b)  the  center  of  pressure 
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pressure  moves  forward  as  the  angle  of  attack  is  increased  (from  0.74c  to  0.31c  for  the  smooth 
airfoil).  As  shown  in  Figure  III.7a,  the  addition  of  the  unactuated  JIF  (which  increases  lift) 
moves  the  center  of  pressure  aft  (0.50c  at  Cl  =  0.4)  and  the  full  range  of  actuation  moves  it 
somewhat  aft  (0.38c  at  Cl  =  0.4)  relative  to  the  smooth  airfoil.  Because  of  the  increase  in  lift  that 
is  associated  with  the  placement  of  the  JIF  on  the  smooth  airfoil,  the  moment  coefficient  is  more 
negative  (compared  to  the  smooth  airfoil)  indicating  increased  stability  (Figure  ni.7a). 
Apparently,  the  suction  that  is  induced  by  the  actuation  jet  near  the  trailing  edge  of  the  airfoil 
contributes  to  a  pitch  up  moment  and  therefore  to  a  reduction  in  Cm  although  it  is  still  negative 
over  the  entire  range  of  angles  of  attack. 

In  general,  use  of  the  actuation  at  the  maximum  attainable  Cft  results  in  a  decrease  in  lift  and  a 
larger  decrease  in  pressure  drag,  resulting  in  an  increased  L/Dp  ratio.  The  effects  of  actuation  are 
also  evident  in  the  pressure  distribution,  which  shows  a  low-pressure  region  upstream  of  the  jet. 
In  addition,  the  increased  suction  and  pressure  that  appear  on  the  suction  and  pressure  sides  of 
the  airfoil,  respectively,  when  the  anactuated  JIF  is  used,  are  diminished  when  actuation  is 
applied. 

The  effect  of  the  streamwise  placement  of  the  JIF  on  the  lower  (pressure)  surface  of  the  airfoil 
was  Investigated  over  a  range  of  angles  of  attack  (for  Cp  =  8x10'^).  Figures  111.8a  and  b  show  the 
variation  of  the  lift  and  pressure  drag  with  (normalized)  xg/c  for  a=Q°  and  8°.  These  data  show 
that  for  xq/c  <  0.7  the  effect  of  the  JIF  (in  the  presence  or  absence  of  actuation)  is  small  at  both 
angles  of  attack.  Compared  to  the  smooth  airfoil,  the  presence  of  the  JIF  within  this  streamwise 
domain  results  in  small  reduction  in  Cl  and  Cop-  Near  the  trailing  edge  of  the  airfoil,  in  the 
absence  of  actuation  the  JIF  leads  to  an  increase  in  both  lift  and  pressure  drag.  The  addition  of 
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Figure  111.8.  The  variation  of  the  (a)  lift  and  (b)  pressure  drag  with  xc/cfor  a=(f  (O)  and  5“  (^).  Solid 
and  open  symbols  correspond  to  the  actuated  and  the  unactuated  cases,  respectively. 
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actuation  results  in  a  decrease  in  lift  (which  is  almost  invariant  with  angle  of  attack)  and  a 
reduction  in  pressure  drag.  However,  the  reduction  in  pressure  drag  as  a  result  of  the  actuation 
increases  with  angle  of  attack  (46%  and  29%  at  or  =  0°  and  8°,  respectively).  It  is  noteworthy  that 
for  xc^c  <  0.75,  the  actuation  results  in  a  slight  increase  in  pressure  drag  ostensibly  as  a  result  of 
the  suction  that  is  induced  by  the  presence  of  the  actuation  jet. 

The  flow  field  in  the  vicinity  of  the  airfoil's  trailing  edge  is  investigated  using  particle  image 
velocimetry  (PIV)  in  the  x-y  plane  (z  =  0).  Each  measurement  field  is  comprised  of  7  partially 


Figure  IIL9.  Global  views  oj 
superimposed  time-averaged  maps  oj 
spanwise  vorticity  and  cross-stream 
velocity  vectors  for  the  smooth  airfoil  (a), 
and  with  the  JIF  in  the  absence  (b)  and 
presence  (c)  of  actuation. 


x/c 


Figure  IILIO.  As  in  Figure  I1L9, 
magnified  views  of  superimposed  time- 
averaged  maps  of  spanwise  vorticity  and 
cross-stream  velocity  vectors. 
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overlapping  frames  measuring  5.56cm  x  5.61cm  where  each  frame  is  computed  from  an 
ensemble  of  500  realization  pairs.  Figures  III.9a-c  show  composite  views  of  time-averaged  maps 
of  cross  stream  vorticity  (that  also  include  regularly-spaced  cross  stream  velocity  vectors 
profiles)  in  the  streamwise  domain  0.95  <x4?  <  1.14  for  the  smooth  (baseline)  airfoil  and  for  the 
airfoil  with  the  mounted  JIF  in  the  absence  and  presence  of  actuation.  In  the  present 
experiments,  the  laser  illumination  is  provided  from  above,  and  therefore  the  flow  below  the 
airfoil  (which  also  includes  the  JIF)  is  shadowed  by  the  airfoil  and  is  not  visible  here.  The 
outline  of  the  airfoil  surface  near  the  trailing  edge  is  included  in  each  frame. 

The  evolution  of  the  near  wake  of  the  smooth  airfoil  is  shown  for  reference  in  (Figure  in.9a). 
The  boundary  layers  on  the  top  and  bottom  surfaces  (marked  by  clockwise  (CW)  and  counter 
clockwise  (CCW)  vorticity  concentrations,  respectively)  merge  downstream  of  the  trailing  edge 
and  commensurate  with  the  (small)  lift  at  a  =  0°,  the  center  of  the  wake  (as  marked  by  local 
velocity  minima)  is  turned  downward.  The  presence  of  the  JIF  (xg/c  =  0.94)  results  in  cross¬ 
stream  spreading  of  the  wake  and  stronger  downward  tilt  (Figure  III.9b).  Compared  to  the 
smooth  airfoil's  wake,  the  (CCW)  vorticity  layer  near  the  lower  (pressure)  surface  of  the  airfoil  is 
thicker  and  appears  to  be  more  intense.  As  shown  in  Figure  III.  11  (below)  the  presence  of  the 
unactuated  JIF  leads  to  organized  (time-periodic)  shedding  of  vortices  that  scale  with  the 
thickness  of  the  (CCW)  vorticity  layer  (nominally  0.02c)  near  the  trailing  edge)  and  apparently 
merge  into  the  airfoil's  wake.  One  of  the  prominent  features  of  the  near  wake  in  Figure  III.9b  is 
the  local  increase  in  its  cross-stream  width  immediately  downstream  of  the  trailing  edge 
compared  to  the  smooth  airfoil.  A  magnified  view  of  the  domain  1  <x4^<  1.03  is  shown  in 


Figure  III.ll.  Power  spectra  of  the  velocity  magnitude  within  the  wake  at  x/c  =  1.03  (a)  and  1 .55  (b)for  the 
smooth  airfoil  (  ),  and  for  the  unactuated  (-  -  -)  and  actuated  JIF  ( - )  configurations. 
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Figure  III.  10b  (the  magnified  view  for  the  smooth  airfoil  is  shown  for  reference  in  Figure 
III.  10a).  This  image  shows  that  the  flux  of  vorticity  from  the  top  and  bottom  eliminated  and  the 
flow  is  turned  strongly  upward. 

Power  spectra  of  the  velocity  magnitude  measured  (using  a  single  hot  wire  sensor)  within  the 
wake  at  x/c=  1.03  and  1.55  are  shown  in  Figures  III.  11a  and  b,  respectively  for  the  smooth 
airfoil,  and  for  the  unactuated  and  actuated  JIF  configurations.  The  spectra  were  measured  on 
both  sides  of  the  wake  at  cross-stream  elevations  where  the  velocity  deficit  is  half  the  maximum 
deficit.  For  clarity,  the  spectra  measured  at  the  lower  part  of  the  wake  are  shifted  vertically 
downwards  by  three  decades.  The  power  spectra  of  the  smooth  airfoil  are  virtually  featureless 
and  are  shown  for  reference  (gray  curve)  in  each  of  the  plots.  The  spectra  of  the  unactuated  JIF 
flow  (dashed  line)  atx/c  =  1.03  (Figure  III.  11  a)  exhibit  a  strong  spectral  peak  at  128  Hz  on  both 
sides  of  the  wake.  This  peak  is  apparently  associated  with  the  shedding  frequency  of  an 
organized  vortex  train  that  is  induced  by  the  JIF  (on  the  pressure  side)  just  upstream  of  the 
trailing  edge.  (The  spectrum  at  the  upper  part  of  the  wake  also  includes  a  weak  spectral  peak  at 
the  first  harmonic)  Similar  vortex  shedding  in  the  presence  of  a  conventional  Gurney  flap  at  the 
trailing  edge  of  an  Eppler  e423  airfoil  was  observed  by  Jeffrey  el  al.'^.  It  is  noteworthy  that 
while  there  is  an  overall  reduction  in  the  magnitude  of  all  spectral  components  at  x/c=  1.55 
(Figure  III.  11b),  the  magnitude  of  the  spectral  peak  at  the  shedding  frequency  remains  relatively 
unchanged  relative  to  the  other  (background)  components. 

The  changes  in  the  power  spectra  in  the  presence  of  actuation  are  remarkable  in  that  the  spectral 
component  at  the  vortex  shedding  frequency  vanishes  completely  and  a  spectral  component  at  the 
actuation  frequency  (1770  Hz)  appears  within  the  dissipation  range  of  the  flow.  At  x/c  =  1.03, 
this  peak  is  somewhat  weaker  on  the  upper  side  of  the  wake,  but  by  x/c  =  1.55  (Figure  III.  1  lb) 
the  magnitude  of  this  peak  is  almost  equal  on  both  sides  of  the  wake  (although  smaller  than  at  the 
low  side  of  the  wake  at  x/c  =  1 .03).  Another  feature  of  the  spectra  of  the  actuated  flow  at 
x/c=  1.03  is  that  the  magnitude  of  all  spectral  components  (and  ostensibly  of  the  total  kinetic 
energy)  is  considerably  smaller  than  the  corresponding  components  in  the  absence  of  actuation 
(i.e.,  with  the  JIF  alone).  In  fact,  with  the  exception  of  the  spectral  peak  at  the  actuation 
frequency,  the  spectra  of  the  smooth  airfoil  and  of  the  actuated  flow  are  virtually  identical  on  the 
top  side  of  the  wake.  However,  on  the  bottom  side  of  the  wake  the  spectral  components  of  the 
actuated  flow  are  larger  than  the  corresponding  components  of  the  smooth  flow  (although  still 
smaller  than  for  the  unactuated  flow).  As  expected,  note  that  the  kinetic  energy  that  is  imparted 
to  the  flow  by  the  actuation  is  primarily  concentrated  within  a  frequency  band  around  the 
actuation  frequency.  The  increase  in  the  magnitude  of  spectral  components  above  the  actuation 
frequency  suggests  enhanced  energy  transfer  at  these  scales  and  therefore  enhanced  dissipation. 
Figure  III.l  lb  shows  probably  the  most  profound  effect  of  the  actuation  on  the  wake  flow.  At 
this  location,  all  the  spectral  components  in  the  actuated  flow  (except  at  the  actuation  frequency) 
are  lower  (on  both  sides  of  the  wake)  than  in  the  wakes  of  the  smooth  airfoil  and  the  airfoil  with 
the  unactuated  JIF.  It  is  remarkable  that  the  effect  of  the  actuation  reduces  the  power  (and 
turbulent  kinetic  energy)  at  both  low  and  high  frequencies  suggesting  that  the  actuation  affect  the 
entire  energy  cascade  within  the  flow  and  effects  enhanced  dissipation. 
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III.6  CONCLUSIONS 

The  objective  of  the  control  approach  described  here  is  to  reduce  the  pressure  drag  at  low  angles 
of  attack  with  minimal  or  no  penalty  in  lift.  The  pressure  distribution  on  the  surface  of  the  airfoil 
is  modified  using  virtual  aero-shaping  that  is  achieved  by  the  formation  of  a  stationary 
recirculating  flow  domain  next  to  the  surface  as  a  result  of  the  interaction  between  the  cross  flow 
and  a  high-frequency  synthetic  jet  actuator  placed  downstream  from  a  miniature  surface-mounted 
passive  obstruction.  The  formation  of  this  interaction  domain  alters  the  flow  above  the  airfoil  by 
inducing  a  displacement  of  the  local  streamlines  that  is  sufficient  to  modify  the  surface  pressure 
distribution  even  at  small  angles  of  attack. 

In  the  present  work,  the  height  of  the  obstruction  (above  the  surface  of  the  airfoil)  is  /i  =  9.3  10'  c 
(i.e.,  less  than  1%  of  the  chord)  and  it  is  placed  at  x</c  =  0.22  where  the  synthetic  jet  actuator  is 
located  immediately  downstream  of  the  obstruction.  The  experiments  are  conducted  at 
Rbc  =  381,000  for  which  the  baseline  flow  is  fully  attached  to  the  surface  of  the  airfoil  for  a <  8“. 

The  interaction  domain  between  the  actuator  and  the  cross  flow  is  investigated  using  high- 
magnification  particle  image  velocimetry  in  the  x-y  plane  (z  =  0).  The  characteristic  streamwise 
scale  of  the  induced  recirculating  flow  domain  is  varied  by  changing  the  actuation  frequency 
is  varied  between  1.25  and  7.8)  and/or  the  jet  strength  (C^,  varies  between  1.7  10*^  and  1.210'0. 
The  present  work  shows  that  the  streamwise  length  of  the  interaction  domain  and  its  cross  stream 
penetration  increase  with  decreasing  actuation  frequency  and  jet  momentum  coefficient.  Phase- 
averaged  PIV  images  show  that  for  =  1.2  10■^  at  given  actuation  frequency  the  streamwise 
length  of  the  interaction  domain  is  typically  comparable  to  one  to  two  wavelengths  of  the 
actuation  frequency.  The  interaction  domain  contains  one  strong  clockwise  vortex  and  a  second 
weaker  vortex  (from  the  earlier  actuation  cycle)  is  advected  downstream  and  rapidly  loses  its 
identity.  However,  a  reduction  in  C^t  by  an  order  of  magnitude  (to  1.710"^)  significantly 
lengthens  the  interaction  domain  indicating  that  it  does  not  merely  scale  with  the  actuation 
frequency.  For  example,  at  =  7.8  (850  Hz)  the  length  of  the  interaction  domain  increases 
from  \.\  Xact  to  iXact  and  at  lower  actuation  frequency  for  which  is  1.25  the  decrease  in  Q 
results  in  a  mismatch  between  the  interaction  domain  and  flow  separation. 

As  a  result  of  the  changes  in  the  curvature  of  the  interaction  domain,  the  induced  low-pressure 
(suction)  peak  on  the  airfoil  diminishes  and  widens  with  the  actuation  frequency  and  jet 
momentum  coefficient  leading  to  a  substantial  increase  in  pressure  drag  and  a  reduction  in  lift.  It 
is  remarkable  that  a  significant  decrease  in  the  baseline  pressure  drag  (compared  to  its 
magnitude)  with  minimal  penalty  in  lift  is  only  achieved  at  higher  actuation  frequencies 
(fact  =  850  Hz).  The  present  results  also  show  that  at  high  actuation  frequencies  the  aerodynamic 
forces  (and  the  shape  of  the  interaction  domain)  are  substantially  less  sensitive  to  changes  in  the 
jet  momentum  coefficient. 
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IV.  INTERACTION  OF  A  SYNTHETIC  JET  WITH  A  TURBULENT  BOUNDARY  LAYER 
IV,1  OVERVIEW 

The  interaction  of  a  synthetic  jet  with  a  cross-flow  fully  developed  turbulent  boundary  layer  over 
a  flat  plate  is  investigated.  Unlike  the  interaction  of  the  jet  with  the  cross  flow  over  an 
aerodynamic  body,  in  the  present  configuration,  the  coupling  between  local  interaction  and 
global  aerodynamic  effects  is  removed.  The  focus  of  this  work  is  the  documentation  of  the  most 
salient  flow  features  of  two  distinct  flow  regimes  where  the  jet  scales  either  with  the  inner 
(viscous)  sub-layer  or  with  the  entire  boundary  layer.  The  streamwise  pressure  gradient,  which 
is  nominally  zero,  is  invariant  throughout  the  experiment. 

IV.1.  FLAT  PLATE  BOUNDARY  LAYER  MODEL 

The  flat  plate  model  forms  one  of  the  short  walls  of  a  25.4  x  45.7  cm  test  section  of  a  closed- 
return  wind  tunnel.  In  order  to  provide  complete  optical  access,  the  working  surface  is  displaced 
into  the  test  section  by  5.1  cm  from  the  original  plane  of  the  tunnel  wall  thus  reducing  the  cross- 
sectional  area  to  25.4  x  40.6  cm.  A  fairing  section  transitions  the  flow  from  outlet  of  the 
contraction,  and  occupies  the  first  26.7  cm  of  the  test  section.  The  flat  plate  model  includes  a 
10.2  cm  long  spanwise  suction  module  immediately  downstream  of  the  leading  edge  fairing 
where  the  flow  rate  (5.2  x  10'^  mVsec  maximum)  is  controlled  by  a  needle  valve  and  monitored 
by  a  Meriam  Laminar  Flow  Element.  The  plate  allows  for  modular  installation  of  the  actuator  in 
order  to  vary  its  streamwise  position  and  is  instrumented  with  static  pressure  ports.  The  two 
synthetic  jet  actuator  modules  were  used.  One  module  is  driven  by  eight  piezoelectric  disks  and 
has  an  orifice  width  of  0.51  mm,  while  the  second  module  is  driven  by  two  loudspeakers  and  has 
an  orifice  width  of  1.91  mm.  Both  jet  orifices  span  the  center  15.2  cm  of  the  test  section. 

In  all  experiments,  the  boundary  layer  was  tripped  by  a  3  cm  long  strip  of  60-grit  sandpaper  that 
was  fastened  to  surface  at  the  contraction  outlet.  The  thickness  of  the  boundary  layer  could  be 
increased  by  a  3.2  mm  tall  fence  that  was  placed  at  the  downstream  edge  of  the  suction  section. 
This  design  enables  variations  of  the  boundary  layer  thickness  5  (at  the  orifice  of  the  jet 
actuator)  between  10  and  40  mm. 

IV.2  CHARACTERIZATION  OF  THE  BASELINE  BOUNDARY  LAYER 

In  order  to  ensure  that  the  boundary  layer  at  the  measurement  station  is  fully  turbulent,  a  strip  of 
60-grit  sandpaper  is  installed  at  the  outlet  of  the  contraction.  Normalized  boundary  layer  profiles 
for  the  baseline  flow  are  shown  in  Figure  IV. la  at  two  different  streamwise  locations  {X=  406 
mm  &  762  mm,  measured  from  the  outlet  of  the  contraction).  There  is  little  variation  in  the 
normalized  profile  for  5  <  C/«,  <  12.5  m/s.  The  boundary  layer  profiles  at  the  conditions  at  which 
synthetic  jet  interaction  is  studied  are  plotted  in  wall  units  (y*,  in  Figure  IV. lb  in  order  to 
compare  to  the  universal  distribution  proposed  by  Coles  (1956)  for  the  velocity  profile  outside  of 
the  viscous  sub- layer:  =\/ Klny^ +C+(n / K)OJ(y/ 5)  .  Where  Lt  =  ujUj  and 

/  =  yUj-fv,  and  the  friction  velocity  (Ut)  is  defined  by  the  wall  shear  stress,  Ut  = 

The  constants  «•=  0.41  and  C  =  5.0,  are  generally  considered  to  be  invariant  with  respect  to 
Reynolds  number  (Erm  &  Joubert,  1991).  The  third  term  in  the  above  equation  is  known  as  “the 
law  of  the  wake”  and  applies  primarily  to  the  outer  region  of  the  boundary  layer.  For  the 
experimental  conditions  used  in  §§IV.3,  IV.4  and  IV.5  {Ree—  520,  2,100,  &  710, 
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Figure  IVJ.  a)  Normalized  velocity  profiles,  (a)  x 
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respectively),  it  is  apparent  that  the  present  velocity  distributions  include  a  logarithmic  region. 
At  Reff=520  and  2,100  reasonable  agreement  with  the  universal  distribution  is  found  for 
*•=  0.41  and  C  =  5.0  with  C/7-=  0.106,  and  0.445  m/s,  respectively.  At  Ree=  710,  the  data  was 
found  to  fit  the  universal  distribution  with  Ut—  0.405  m/s  with  K=  0.41  and  C  =  4.7.  Within  the 
viscous  sub-layer  if  where  the  edge  of  this  part  of  the  boundary  layer  occurs  near  3;^  ~  15. 
Cross-stream  distributions  of  Reynolds  stresses  within  the  baseline  boundary  layers  are  shown  in 
Figure  1V.2.  The  distributions  at  the  three  Reynolds  numbers  are  qualitatively  similar  and 
magnitude  of  the  peak  stresses,  which  occur  near  the  inner  edge  of  the  logarithmic  region,  are  in 
agreement  with  the  levels  within  fully  turbulent  boundary  layers  at  similar  Reynolds  numbers 
reported  by  Erm  &  Joubert  (1991).  Based  on  these  comparisons  of  the  mean  streamwise  velocity 
and  turbulent  stresses,  the  baseline  boundary  layers  in  the  present  investigation  are  judged  to  be 
fully  turbulent.  A  detailed  discussion  of  the  characteristics  of  low  Reynolds  number  {Ree  < 
6,000)  turbulent  boundary  layers  downstream  from  various  tripping  devices  can  be  found  in  a 
paper  by  Erm  &  Joubert  (1991),  that  also  includes  an  extensive  review  of  earlier  work. 

IV.3  CLOSED  RECIRCULATION  REGIONS 

As  shown  in  earlier  work,  the  interaction  between  a  synthetic  jet  and  cross-flow  leads  to  the 
formation  of  a  closed  recirculation  region  if  the  dimensionless  formation  frequency,  /,  is  above 
some  critical  value.  This  dimensionless  frequency  is  defined  as  the  product  of  Strouhal  number 
(based  on  jet  orifice  width)  and  dynamic  pressure  ratio  (jet  to  cross-flow).  The  characteristic 
features  of  closed  recirculation  regions  are  investigated  in  turbulent  boundary  layer  where  Ree  = 
500,  &b  =  73,  Sts  =  21,  for  UjIU^  =  3.0,  4.0,  5.0  and  6.0.  (The  dimensional  parameters  are 

Uoo  =  2.0  m/s,  <J=37.0  mm,  ^=4.0  mm,  6  =  0.51  mm, /=  1130  Hz,  and  Uj  =  6.0,  8.0,  10.0  & 

12.0  m/s).  The  local  dimensionless  frequencies  are  /  =  0.6,  1.2, 1.8,  and  2.6,  all  well  above  the 
critical  value  (0.1)  at  which  a  closed  recirculation  region  forms. 

The  time-averaged  flow  field  in  the  vicinity  of  the  synthetic  jet  actuator  orifice  is  measured 
within  twelve  partially  overlapping  sets  of  PIV  images  acquired  at  33.1  jim/pixel  (see  table  A.2). 
Concentration  maps  of  normalized  time-averaged  spanwise  vorticity,  s  are  shown 

in  Figure  IV.3  (Left)  (positive  vorticity  is  CCW).  The  closed  recirculation  regions  are  outlined 
using  the  streamline  ^  =  0,  where  the  streamfunction  is  computed  by  cross-stream  integration  of 
the  streamwise  velocity  component,  with  ^  =  0  at  the  surface  of  the  plate.  It  is  apparent  that  the 
penetration  of  the  synthetic  jet  increases  with  dynamic  pressure  ratio  in  a  similar  manner  to  a 
conventional  jet  in  cross  flow.  The  synthetic  jet’s  counter  rotating  vortex  pairs  appear  as  bands 
of  positive  and  negative  vorticity  that  are  deflected  downstream  by  the  cross  flow  and  dissipate 
rapidly.  A  closed  recirculation  region  forms  as  the  cross  flow  separates  just  upstream  of  the 
actuator  orifice  and  reattaches  farther  downstream  (x/8  =  Ltubbie)  forming  a  new  dividing 
streamline  {y/  =  0)  between  upstream  and  downstream  stagnation  points.  The  characteristic 
scale  of  the  closed  recirculation  increases  with  dynamic  pressure  ratio  {Lb^iJS  =  0.60,  0.95, 1.1 
and  1.2  for  VR  =  3.0,  4.0,  5.0  and  6.0,  respectively),  and  a  similar  dependence  of  has  also  been 
observed  by  Mittal  &  Rampunggoon  (2002).  During  each  suction  cycle  the  synthetic  jet  actuator 
draws  in  fluid  originating  both  upstream  and  downstream  from  the  orifice,  removing  some  fluid 
from  the  recirculation  bubble,  and  subsequently  replacing  it  during  the  blowing  cycle.  A  simple 
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analysts  shows  that  the  total  volume  of  fluid  moved  through  the  actuator  orifice  during  each  half¬ 
cycle  is  less  than  2  %  of  the  volume  of  the  recirculation  region  at  every  velocity  ratio.  Thus,  the 
size  of  the  recirculation  region  does  not  vary  significantly  over  the  actuation  cycle.  As  the 
velocity  ratio  is  increased,  the  aspect  ratio  of  the  recirculation  region  (length-to-height)  changes 
from  approximately  4  (UjfU^  =  3.0  &  4.0)  to  2  (Uj  =  5.0  and  6.0). 

The  presence  of  the  recirculation  region  clearly  distorts  the  velocity  field  and  alters  the 
downstream  evolution  of  the  boundary  layer.  Distributions  of  the  mean  streamwise  velocity  are 
shown  in  Figure  IV.3-Right  at  three  streamwise  locations  for  each  velocity  ratio.  Upstream  from 
the  actuator  orifice  (x/S  =  -0.5)  the  boundary  layer  is  thicker  than  in  the  baseline,  due  to  the 
blockage  that  is  introduced  by  the  synthetic  jet.  At  ji:/^=  0.5,  there  is  a  local  velocity  maximum 
near  the  outer  edge  of  the  recirculation  region  (defined  by  if^=  0),  which  increases  with  velocity 
ratio.  At  higher  cross-stream  elevations  the  streamwise  velocity  is  higher  than  the  baseline  at  all 
velocity  ratios,  indicating  that  the  cross  flow  is  accelerated  due  to  the  displacement  of  the 
recirculation  region.  At  this  streamwise  location  the  strong  reverse  flow  with  the  recirculation 
region  toward  the  actuator  orifice  is  also  evident.  While  downstream  from  the  interaction 
domain  the  CW  vorticity  is  quickly  reorganized  into  a  new  inner  layer  immediately  adjacent  to 
the  surface  [cf.  Figure  IV.3  (Left)],  the  return  to  the  velocity  distribution  of  the  baseline  flow 
requires  a  much  greater  streamwise  distance.  At  x/S=  3,  an  inner  region  of  the  boundary  layer  is 
re-established  for  all  velocity  ratios,  however  the  velocity  distributions  are  still  substantially 
distorted  relative  to  the  baseline  flow.  This  is  particularly  evident  as  the  velocity  ratio  is 
increased,  due  to  the  wake  effect  (or  velocity  deficit)  that  is  associated  with  the  blockage  caused 
by  the  recirculation  region.  The  velocity  distribution  for  VR  =  6.0  has  a  local  maximum  at  the 
edge  of  the  new  inner  layer,  followed  by  a  velocity  deficit  at  higher  elevations  before  reaching 
the  free-stream  velocity.  It  is  noteworthy  to  mention  that  CCW  vorticity  is  generated  in  the 
region  between  the  outer  jet  flow  and  the  newly  formed  inner  layer  downstream  from  the 
reattachment  point  of  the  recirculation  region  [Figure  IV.3  (Left)].  The  rapid  development  of  the 
inner  layer  is  also  evident  in  the  normalized  wall  shear  stress,  =  T„l{^pUi),  shown  in  Figure 
IV.4-Left.  At  every  velocity  ratio,  the  blockage  caused  by  the  recirculation  region  causes  lower 
shear  stress  upstream  from  the  actuator  orifice,  compared  to  the  baseline  value.  The  shear  stress 
is  negative  within  the  recirculation  region,  and  then  rapidly  rises  as  the  boundary  layer 
redevelops  downstream.  At  every  velocity  ratio,  the  shear  stress  reaches  a  maximum  value 
within  1.5  ^from  the  reattachment  point,  confirming  the  rapid  formation  of  the  inner  wall  region 
within  the  boundary  layer.  Farther  downstream,  the  streamwise  variation  of  the  shear  stress  is 
reduced  as  the  velocity  distribution  relaxes  to  that  of  the  baseline  boundary  layer. 

The  effect  of  the  presence  of  the  closed  recirculation  regions  on  the  boundary  layer’s 
displacement  thickness  momentum  thickness  6,  and  shape  factor  1^2  =  S'  are  shown  in 
Figure  IV.4-Right.  As  the  velocity  ratio  is  increased,  the  displacement  thickness  within  the 
interaction  domain  increases  due  to  the  increased  height  of  the  recirculation  region  [cf.  Figure 
IV.3  (Left)].  Concomitantly,  the  convective  acceleration  of  cross-flow  fluid  and  the  additional 
momentum  injected  by  the  synthetic  jet  result  in  substantial  reduction  of  the  momentum 
thickness  within  the  interaction  domain.  The  momentum  thickness  becomes  negative  for 
UjU„=  4.0,  5.0  &  6.0,  indicating  that  the  flow  contains  more  momentum  than  a  uniform 

profile  at  f/«.  Downstream  from  the  recirculation  region,  the  increased  momentum  thickness 
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compared  to  the  baseline  flow  is  attributed  to  the  wake  effect  that  is  evident  in  the  velocity 
distribution  (Figure  IV.3-Right).  The  shape  factor  is  reduced  since  a  significant  amount  of 
momentum  is  transported  into  the  newly  formed  inner  layer  (Figure  IV.3-Right). 

IV.4  JET  PENETRATION  SCALING  WITH  THE  INNER  LAYER 

In  this  section,  the  interaction  of  a  synthetic  jet  with  the  flat  plate  boundary  layer  is  studied  under 
conditions  leading  to  the  formation  of  discrete  vortices  that  scale  with  the  inner  segment  of  a 
turbulent  boundary  layer.  This  investigation  is  carried  out  at  Reg  =  2,100,  ^  =  61, 
UjIu^  =  1.25,  and  Sts  =  3.5  (f/„  =  10  m/s,  S=  31.0  mm,  9=  3.2  mm,/=  1130  Hz,  U 12.5 

m/s,  and  b  =  0.051  mm).  The  local  dimensionless  frequency  at  these  parameters  is  /=0.02 
(assuming  the  mean  jet  centerline  velocity  in  the  absence  of  cross-flow  106  from  the  exit  plane  is 
--  0.5  Uj ).  Phase-locked  measurements  are  obtained  at  eighteen  instances  during  the  actuation 

cycle  (in  20°  increments)  using  ten  overlapping  sets  of  PIV  images.  Six  sets  of  400  images  per 
phase  were  acquired  at  10.4  |im/pixel  to  resolve  the  flow  immediately  adjacent  to  the  wall,  and 
four  sets  of  100  images  per  phase  were  acquired  at  a  resolution  of  46.9  ixm/pixel  to  resolve  the 
outer  region  of  the  boundary  layer. 

IV.4.1  The  Velocity  Field 

The  interaction  between  the  synthetic  jet  and  cross  flow  leads  to  the  formation  of  discrete 
vortices  having  a  characteristic  cross  section  that  is  smaller  than  the  boundary  layer  thickness 
(0.08  <5),  as  shown  in  the  phase-averaged  maps  of  spanwise  vorticity  concentrations  in 
Figure  IV.5.  The  origin  of  the  coordinate  system  is  at  the  centerline  of  the  actuator  exit  plane 
and  ^  =  0°  corresponds  with  the  beginning  of  the  blowing  phase  of  the  cycle.  Selected  cross¬ 
stream  velocity  distributions  are  also  included.  The  formation  of  the  vortex  pair  during  the 
blowing  phase  of  the  cycle  is  visible  during  0°  <  (p  <  180°  (Figure  IV.5a-c).  A  fraction  in  the 
fluid  expelled  from  the  actuator  rolls  up  into  the  vortex  pair,  resulting  in  a  stem  of  jet  fluid 
behind  the  vortex.  The  upstream  (CCW)  vortex,  which  is  advected  by  faster  moving  flow, 
overtakes  the  CW  vortex.  During  the  suction  phase  (180°  <  (p  <  360°,  Fi^re  IV.5d-f)  the 
boundary  layer  upstream  from  the  actuator  orifice  becomes  thinner  as  fluid  is  drawn  into  the 
actuator.  The  vortices  quickly  diminish  in  strength  and  completely  disappear  within  three 
actuation  wavelengths,  indicating  that  there  is  a  finite  domain  where  the  vortices  directly 
influence  the  flow.  In  the  time-averaged  flow,  the  finite  interaction  domain  appears  as  a 
stationary  vortex  adjacent  to  the  surface  and  immediately  downstream  from  the  actuator  orifice, 
as  shown  in  Figure  IV.6b  (the  baseline  flow  shown  in  Figure  IV.6a).  The  time-averaged  flow 
bears  a  remarkable  resemblance  to  the  quasi-steady  closed  recirculation  regions  discussed  in  the 
previous  section.  The  height  and  width  of  the  interaction  domain  (as  characterized  by  the 
vorticity  concentration)  are  approximately  0.08^  and  0.3^,  respectively.  Downstream  from  the 
interaction  domain,  beginning  eXx/d=  0.4,  the  vorticity  is  reorganized  into  a  thin  layer  adjacent 
to  the  surface. 

Cross-stream  distributions  of  the  time-averaged  streamwise  velocity  at  four  streamwise  locations 
{x/8=  0.13,  0.25,  0.5  &  1.0)  are  shown  in  wall  units  in  Figure  IV.7  using  the  friction  velocity  of 
the  baseline  flow  {Uj  =  0.445  m/s).  Immediately  downstream  from  the  actuator  orifice  (x/^= 
0.13),  the  passage  of  discrete  vortices  results  in  a  local  maximum  aty"^  =  70,  and  negative 
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Figure  IV.3.  Left:  Normalized  spanwise  vorticity,  and  velocity,  ReO-  520,  StS  =  21,  3'b  =  73:  (a)  baseline,  (b)  V jjU,,  ~ 
(c)  4.0,  (d)  5.0,  and  (e)  6.0.  y/=0  is  superimposed,  and  Right:  Mean  velocity:  baseline  ( O),  UjfU„  =  3.0  (0),  4.0  (A),  5.0 
(K),  and 6.0  (T).  x/S  =-0.5,  0. 5,  and 3. 0  (from  top  to  bottom). 


x/5 


FigureIV.4.  Left:  Normalized  wall  shear  stress,  UjlU,,  =3.0(0),  4.0(A),  5.0  (M),  and6.0(r).  The  baseline  stress^  is 
indicated  using  a  dashed  line.Right:  Boundary  layer  displacement  thickness,  momentum  thickness  and  shape  factor:  Uj/U„ 
3.0  (0),  4.0  (A),  5.0  (m,  and  6.0  (n  The  baseline  quantities  are  shown  by  the  dashed  lines. 


velocity  adjacent  to  the  wall.  Beyond  y*  =  200  there  is  little  influence  on  the  mean  velocity 
profile.  At  x/S=  0.25  the  mean  velocity  adjacent  to  the  wall  is  significantly  lower  than  in  the 
baseline,  but  not  reversed.  The  displacement  of  the  velocity  maximum  induced  by  the  vortices  to 
a  slightly  higher  elevation  =  80)  indicates  the  cross  stream  spreading  of  their  cores  as  they 
travel  downstream.  This  trend  continues  at  x/S=  0.5,  with  a  slight  perturbation  associated  with 
velocity  induced  by  the  CW  vortices  located  at  =  100.  In  the  region  dovmstream  from  the 
interaction  domain,  the  velocity  distributions  are  multi-tiered,  with  an  inner  region  adjacent  to 
the  wall  (y^  <  20,  at  x/S=  0.5),  a  middle  region  that  contains  the  remnants  of  vorticity  from  the 
coherent  structures  (20  <  200),  and  finally  an  outer  region  of  the  boundary  layer  (y"'  >  200), 

which  appears  to  be  unaffected  by  the  interaction. 

Since  the  vortices  formed  by  the  interaction  rapidly  diminish  in  strength  as  they  are  advected 
downstream,  the  largest  velocity  fluctuations  are  confined  to  the  interaction  domain,  as 

illustrated  in  the  three  normalized  components  of  Reynolds  stress  {u'u'/ul, ,  v'v'fui  ,  and 

uVlul )  shown  in  Figure  IV.8.  Within  this  domain,  all  three  components  of  Reynolds  stress  are 
much  higher  (at  least  one  order  of  magnitude)  than  the  corresponding  levels  within  the  baseline 
turbulent  boundary  layer.  Cross-stream  distributions  of  the  normal  and  shear  components  of 
Reynolds  stress  are  shown  at  five  streamwise  locations  {x/S=  0.13,  0.25,  0.5,  1.0,  &  1.5)  in 
Figures  IV.9-Lefl:  and  IV.9-Right,  respectively.  From  this  data  it  is  immediately  apparent  that 
the  Reynolds  stresses  are  rapidly  attenuated  with  downstream  distance  from  the  actuator  orifice, 
returning  to  values  that  are  close  to  those  in  the  baseline  boundary  layer  by  x/S=  1.5. 
Additionally,  the  velocity  fluctuations  induced  by  the  interaction  are  confined  to  the  inner  part  of 

the  boundary  layer  (y/S<  0.4).  The  distribution  of  uV/ul  at  x/S=  0.13  (Figure  IV.9-Righta) 
changes  sign  close  to  wall,  ostensibly  the  result  of  mean  reverse  flow  immediately  downstream 
from  the  actuator  orifice  (cf.  Figure  IV.7). 

IV.4.2  The  Pressure  Field 

The  phase-averaged  pressure  gradients  within  the  interaction  domain  are  computed  using  the 
Navier-Stokes  equations  along  with  triple  decomposition  of  the  velocity  field  data.  The 
normalized,  phase-averaged  cross-stream  ^dp/dy>=S/\pU^<dpldy>\  and  streamwise 
<3/7/3x>]  gradients  are  shown  for  six  phases  of  actuation  in  Figure 

IV.  10.  Each  discrete  vortex  produces  large  gradients  in  both  coordinate  directions,  which 
diminish  rapidly  with  distance  from  the  actuator  orifice,  further  indicating  the  finite  nature  of  the 
interaction  domain.  The  time-averaged  modification  of  the  pressure  gradients  within  this  region 
are  shown  in  FigureIV.il,  computed  from  the  Reynolds-averaged  Navier-Stokes  equations. 
Notice  that  the  time-averaged  gradients  within  the  interaction  domain  (0  <  x/S  <  0.3)  have  a 
qualitatively  similar  structure  compared  to  the  phase-averaged  gradients  associated  with  a  single 
vortex.  Thus,  although  the  time-averaged  flow  does  not  exist  in  reality,  the  interaction  between 
the  synthetic  jet  and  cross  flow  leads  to  the  formation  of  a  time-averaged  stationary  vortex.  The 
time-averaged  normalized  pressure,  plotted  in  Figure  IV.  12,  is  computed  using  a  combined-path 
Runge-Kutta  integration  scheme,  with  the  boundary  condition  p(x,<»)  =  0.  (The  pressure  is 
computed  as  gage  pressure  using  the  local  static  pressure  as  the  reference.)  The  lowest  pressure, 
p  =  -0.55(^  pUi ),  in  the  flow  is  near  the  center  of  the  time-averaged  stationary  vortex  {x/S= 
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Figure  IV.  7.  Time-averaged  velocity  profiles  in  wall  units.  From  top  to  bottom,  x/S  -0.13,  0.25,  0.5,  &  1.0, 
with  successive  locations  separated  by  one  decade. 
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Figure  IV.8.  nujul, ,  vv'lul^  and  uvjul,  baseline  (a),  (b),  and  (c),  and  actuated  (d),  (e)  and  (f). 
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Figure  IV.IO.  Phase-averaged  normalized  pressure  gradients,  <dpldx>  (a)-(f);  <dpldy>  (g)-(l):  (a,  g)  0°, 
(b,  h)  60  °  (c,  i)  120  ‘’.{cl,  1)180°,  (e,  k)  240  °,and  (f,  1)  300  ® 
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Figure  IV.IL  Normalized  Reynolds-averaged  pressure  gradient  dpjdx  and  dpl^y  baseline  (a)  &  (b); 
actuated  (c)  &  (d). 
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Figure  IV.  1 2,  Normalized  pressure  field  with  ^  =  0.01  superimposed:  (a)  baseline,  (b)  actuated.  Dashed 
line  indicates  least  significant  contour  of  negative  pressure,  p  =  -0.  015{\pUi) 
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0.\2,  y/S=  0.04).  The  displacement  of  the  cross  flow  around  the  interaction  domain  is  evidenced 
by  the  streamline  y/  =  y//  U^S  =  0.01 . 


IV.43  Boundary  Layer  Evolution 

The  normalized  surface  shear  stress  \pUl)(du /dy)\,  plotted  in  Figure  IV.13,  further 

demonstrates  the  fact  that  the  direct  influence  of  the  discrete  vortices  is  confined  to  a  finite 
domain.  Within  the  interaction  domain,  there  are  large  fluctuations  in  associated  with  the 
passage  of  coherent  vortices.  The  actuation  wavelength  (distance  between  successive  vortices)  is 
approximately  A„ct  0A2S,  easily  determined  from  successive  local  minima  in  shear,  indicates 
that  the  vortex  celerity  {Uc  =f^ct)  is  approximately  0.4  Uoc  Since  the  CW  vortices  0)  are 
stronger  and  closer  to  the  surface  (cf.  Figure  IV.5),  there  is  a  reduction  in  the  time-averaged 
shear  stress  within  the  interaction  domain  (0<x/^<0.3).  Further  downstream,  the  time- 
averaged  shear  is  higher  than  that  present  in  the  baseline  due  to  the  formation  of  an  inner  layer  in 
flow  (cf.  Figure  IV.7).  Interestingly,  downstream  from  the  interaction  domain  there  are  small 
phase-averaged  oscillations  in  the  shear  stress  about  the  time-averaged  value,  that  may  be  the 
result  of  pressure  waves.  Notice  that  for  ^=0°,  60°,  and  120°  the  phase-averaged  wall  shear 
stress  is  uniformly  higher  than  the  time-average  (by  about  5  %)  within  the  domain 
(0  <x^^<  0.3).  Within  this  domain,  the  time-  and  phase-averaged  shear  stress  are  approximately 
equal  at  180°,  while  the  phase-averaged  stress  is  uniformly  lower  for  240°,  and  300. 

The  alteration  of  the  boundary  layer  characteristics  due  to  jet  actuation  is  shown  in  Figure  IV.14. 
In  the  baseline  flow,  the  boundary  layer  displacement  and  momentum  thickness  monotonically 
increase  with  downstream  distance.  At  the  location  of  the  actuator  orifice  the  thickness  based  on 
0.99  is  31.0  mm.  The  shape  factor  remains  constant  over  the  measurement  domain  at  hi2  = 
1.38,  a  typical  value  for  a  fully  turbulent  boundary  layer  (Schlichting,  p.  454,  1979).  As  seen  in 
the  velocity  profiles  (cf.  Figure  IV.7),  the  synthetic  jet  produces  large  changes  in  the  boundary 
layer  within  the  interaction  domain  (0  <  x/S  <  0.5).  In  this  region,  the  cross-flow  is  displaced 
(increased  and  accelerated  (reduced  6),  which  brings  about  large  changes  in  the  shape  factor. 
Downstream  from  the  interaction  domain  <5^,  0,  and  hi2  return  to  their  baseline  values  within  one 
boundary  layer  thickness. 


IV.4.4  Turbulent  Kinetic  Energy 


High  frequency  actuation  has  previously  been  shown  to  enhance  dissipation  in  free  shear  flows 
(Wiltse  &  Glezer,  1998).  In  the  bluff  body  aerodynamics  experiments,  some  of  the  reduction  in 
turbulent  kinetic  energy  within  the  separated  shear  layer  may  be  due  to  this  effect.  The  present 
data  allows  direct  computation  of  the  turbulent  kinetic  energy,  as  well  as  the  rates  of  turbulent 
production  and  turbulent  viscous  dissipation  within  the  interaction  domain.  The  governing 
equation  for  turbulent  kinetic  energy  equation  in  Cartesian  coordinates  is: 
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Figure  IVJ3.  Normalized  wall  shear  stress:  (A)  phase-averaged,  from  top  to  bottom,  60 , 120 , 180  , 

240  ^  and  300  —  baseline;  actuated.  Successive  phases  separated  by  0. 01 


where  ={uY  ^  repeated  indices  indicate  summation,  and  over-bars  indicate 

time-averages.  A  derivation  of  this  equation  can  be  found  in  Viscous  Fluid  Flow  by  White 
(1974).  The  left  hand  side  of  the  equation  is  the  total  rate  of  change  of  the  turbulent  kinetic 
energy,  while  the  terms  on  the  right  hand  side  are  the  rates  of  convective  diffusion,  production, 
viscous  work,  and  viscous  dissipation,  respectively.  The  time  averaged  turbulent  kinetic  energy 
q,  rate  of  turbulence  production  P,  and  rate  of  viscous  dissipation  e  are  shown  in  Figure  IV.IS 
for  the  baseline  and  actuated  flows.  These  quantities  are  computed  from  the  PIV  data  assuming 
the  flow  to  be  two-dimensional  (i.e.  wsO,  anda(  )/dz  =  0),  and  are  normalized  such  that 

q^  =r/Ul,  P^-pi^j(duj/dx,)/apUl{Ujd),zx^^ 

e  =  p  [duydxj  +du'j/dx,ydUj/dXi)/(-^pUl{U^fS))  .  It  is  likely  that  the  turbulence 

production  can  be  computed  to  reasonable  accuracy  using  the  above  assumptions,  since  the 
neglected  terms  contain  either  the  mean  spanwise  velocity  component  (  w  ),  or  spanwise 
derivatives  of  the  mean  streamwise  and  cross-stream  velocities  all  of  which  are  likely  to  be  very 
nearly  zero.  Although  w'  is  definitely  non  zero,  and  its  omission  will  likely  lead  to  errors  in  the 
turbulent  kinetic  energy  and  dissipation,  its  contribution  was  not  included. 

In  the  baseline  flow,  turbulent  kinetic  energy  is  distributed  throughout  the  boundary  layer 
(profiles  of  mV  and  77  are  shown  in  Figure  IV.2)  while  the  highest  levels  of  production  and 
viscous  dissipation  are  adjacent  to  the  wall.  In  the  actuated  flow,  the  levels  of  q^ ,  P,  and  i 
within  the  interaction  domain  are  an  order  of  magnitude  greater  than  those  present  in  the  baseline 
flow.  Negative  turbulence  production,  occurring  near  the  actuator  orifice,  is  interpreted  as  a 
transfer  of  energy  from  the  periodic  jet  motion,  into  the  mean  flow.  The  total  amount  of 

turbulent  kinetic  energy  within  the  boundary  layer  ,  as  well  as  the  total  rates  of  P ,  and  f  as  a 
function  of  x  are  determined  by  integrating  these  quantities  in  the  cross-stream  direction. 
Namely, 

1  ®2i(y  ^  1  o.2i<y^  ^  1  o.2i<y 

Q^{x)  =  -  =  ^  \p{x,y)dy,md  E{x)  =  -  \e{x,y)dy. 

These  quantities  are  integrated  across  the  full  extent  of  the  high-resolution  data  (0  ^y/d'^  0.21) 
and  their  streamwise  variation  is  shown  in  Figure  IV.  16  for  the  baseline  and  actuated  flows.  In 

the  baseline  flow,  is  invariant  with  streamwise  distance,  since  the  rates  of  turbulence 
production  and  viscous  dissipation  are  balanced.  The  large  increase  in  turbulent  kinetic  energy 
within  the  interaction  domain  is  immediately  apparent,  and  its  peak  value  occurs  at  x/5=  0.11. 
The  turbulence  production  and  viscous  dissipation  are  also  greatly  elevated  within  the  interaction 
domain.  Downstream  from  the  interaction  domain,  all  three  quantities  rapidly  decay,  and  by 
x/5=  1  have  almost  recovered  their  baseline  values.  However,  in  contrast  to  high-frequency 
excitation  in  free  shear  flows,  the  total  amount  of  turbulent  kinetic  energy  within  the  flow  does 
not  fall  below  the  baseline  value. 
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The  rapid  attenuation  of  turbulent  kinetic  energy  downstream  from  the  interaction  domain  is  also 
apparent  in  power  spectra  of  the  streamwise  velocity  component,  acquired  using  LDV.  The 
power  spectral  density  is  plotted  against  normalized  wavenumber  {k  =  27tf /Uj)  at  constant 

elevation  within  the  boundary  layer  0/^=  0.065)  for  several  streamwise  locations  in 
Figure  IV.  17-Left,  for  the  baseline  and  actuated  flows.  The  slight  upturn  in  the  power  spectra  at 
the  upper  end  of  the  frequency  range  is  due  to  aliasing  of  energy  contained  above  the  Nyquist 
frequency.  Immediately  downstream  from  the  orifice  (x/<?=0.16),  the  spectral  peaks  at  the 
actuation  frequency  (1130  Hz,  kd=^  22.0)  and  its  higher  harmonics  are  present,  and  there  is 
substantially  more  energy  at  all  frequencies  compared  to  the  baseline.  The  inertial  subrange  is 
shifted  to  higher  frequencies,  starting  at  the  actuation  frequency,  as  energy  is  transferred  to  scales 
above  and  below  the  actuation  frequency.  The  rapid  reduction  in  turbulent  kinetic  energy  is 
evident  in  the  spectra  downstream  from  the  interaction  domain  {x/5=  0.65,  1.3, 1.9,  2.6  and  3.5). 
The  spectral  peak  at  the  actuation  frequency  decays  with  downstream  distance  from  the  actuator 
orifice,  and  the  normal  energy  cascade  is  re-established  as  the  spectra  relax  toward  its  baseline 
shape.  It  should  be  noted  that  atx/S=  1.3, 1.9  &  2.6,  there  is  a  range  of  frequencies  centered  on 
kS~  2,  where  the  actuated  spectra  falls  below  the  baseline,  ostensibly  the  result  of  the  enhanced 
dissipation.  Further  downstream  there  is  little  difference  between  the  baseline  and  actuated 
spectra,  except  for  a  small  peak  at  the  actuation  frequency.  The  variation  in  power  spectra  of  u 
across  the  boundary  layer  is  shown  atx/S=  0.65, 1.3,  and  1.9  in  Figures  IV.  17-Right,  IV.  18-Left, 
and  IV.  18-Right,  respectively.  Similar  to  the  observed  changes  in  the  velocity  distributions,  the 
effect  of  synthetic  jet  actuation  on  the  power  spectra  is  confined  to  the  inner  portion  of  the 
boundary  layer  (y/S<  0.2). 

IV.5  JET  PENETRATION  SCALING  WITH  THE  ENTIRE  BOUNDARY  LAYER 

The  interaction  of  a  synthetic  jet  with  a  flat  plate  boundary  layer  is  also  investigated  under 
circumstances  that  lead  to  the  formation  of  discrete  vortices  that  scale  with  the  entire  width  of  the 
boundary  layer.  This  is  not  a  trivial  undertaking,  since  even  with  maximum  suction  applied  at 
the  leading  edge  of  the  flat  plate  model,  the  boundary  layer  thickness  within  the  test  section  is 
greater  than  10  mm.  In  order  to  create  the  largest  possible  discrete  vortices,  the  loudspeaker 
synthetic  jet  actuator  was  employed  since  the  volume  of  fluid  expelled  per  cycle  is  inversely 
proportional  to  formation  frequency.  Several  combinations  of  actuation  frequency,  amplitude, 
and  free  stream  velocity  were  tried  before  settling  on  the  conditions:  t/«=  8.0  m/s,  S=  12.5  mm, 
0=  1.4  mm,  /=  300  Hz,  Uj  =  12.5  m/s,  and  6  =  1.91  mm.  The  corresponding  dimensionless 

parameters  are:  Ree—  700,  &b  —  6.5,  U jJU^  =  1.25,  and  Sts—  0.47.  In  the  context  of  the  flow 
regime  transition  map  (cf.  Figure  22)  the  local  dimensionless  frequency  for  this  set  of  parameters 
is  /  =  0.03  (assuming  the  mean  centerline  velocity  at_y  =  106  is  -0.5  Uj).  Eighteen  phases  of 

the  flow,  locked  to  the  actuation  in  20°  increments,  are  measured  in  fifteen  overlapping  sets  of 
PIV  images.  Ten  sets  of  400-800  images  per  phase  are  acquired  at  15.1  pm/pixel  to  resolve  the 
flow  immediately  adjacent  to  the  wall,  and  five  sets  of  400  images  per  phase  are  acquired  at  30.7 
pm/pixel  resolve  the  outer  region  of  the  boundary  layer. 
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Figure  IVJ7  u power  spectra  (baseline  is  in  gray).  Left:  y/8  =  0. 065  (y-^  =  60),  x/S  =0.16,  0. 65,  1.3,  1.9,  2.6,  and 
3.5  (top  to  bottom).  Right:  x/5=  0.65,  y/S=  0.032,  0.065,  0.097,  and  0.16  (top  to  bottom)  (y+  =  30,  60,  90,  and 
145).  Successive  locations  displaced  by  two  decades. 


Figure  IV,  18.  u power  spectra  (baseline  in  gray)  Left:  x/S=  1.3,  y/8=  0.032,  0.065,  0.097,  0.16,  and  0.26  (top  to 
bottom).  (>+  =  30,  60,  90,  145,  and 240)  Right:  x/8=  1.9,  y/8  =  0.032,  0.065,  0.13,  0.26,  and  0.52  (top  to  bottom). 
(y+  =  30,  60,  120,  240,  and  470.  Successive  locations  displaced  by  two  decades. 


IV.5.1  The  Velocity  Field 

The  phase-averaged  spanwise  vorticity  field  is  shown  for  every  third  phase  of  actuation  in 
Figure  IV.  19-Left.  At  the  onset  of  the  blowing  portion  of  the  cycle  (^=0°,  Figure  IV.  19-Left  a), 
the  previously  formed  vortices  from  the  interaction  have  been  advected  downstream  from  the 
actuator  orifice  by  the  free  stream.  During  the  blowing  cycle  a  vortex  pair  forms  at  the  actuator 
orifice.  By  ^  =  60°  (Figure  IV.  19-Left  b)  the  formation  of  the  vortex  pair  is  almost  complete, 
and  as  the  blowing  cycle  continues  a  stem  of  fluid  is  expelled  behind  the  pair  (Figure  IV.  19-Left 
c).  The  vortex  pairs  penetrate  to  approximately  y/S  =  0.9  at  ^  =  240°,  much  farther  than  the 
previous  case  (§IV.4,  Ree  =  2,100,  S^b  =  61.  =  1.25,  &  Sts=  3.5).  As  the  jet  fluid  is 

deflected  by  the  cross-flow,  the  upstream  (CCW)  half  of  the  vortex  pair  is  subjected  to  a  higher 
advection  velocity  and  overtakes  the  CW  vortex  (Figure  IV.  19-Left  d  and  e).  During  the  suction 
phase  (180°  <360°,  Figure  IV.  19-Left  d-f),  upstream  fluid  is  drawn  into  the  actuator, 

resulting  in  a  “fuller”  boundary  layer  profile  immediately  upstream  from  the  actuator  orifice. 
Once  the  discrete  vortices  are  fully  formed,  they  remain  at  approximately  the  same  elevation 
within  the  boundary  layer  and  rapidly  diminish  in  strength  and  coherence  as  they  travel 
downstream.  The  vortices  completely  disappear  within  three  wavelengths,  and  a  new  inner 
boundary  layer  is  re-established. 

In  the  time-averaged  interaction  domain  is  comprised  of  a  stationary  vortex  adjacent  to  the 
surface  and  immediately  downstream  from  the  actuator  orifice  as  shown  in  the  time-averaged 
vorticity  field  in  Figure  132.  The  vortex  extends  approximately  O.S^outward  from  the  surface 
and  2S  downstream.  The  discrete  vortices  are  of  sufficient  strength  to  induce  time-averaged 
reverse  flow  within  the  interaction  domain,  and  therefore  a  time-averaged  separation  bubble  is 
created  as  shown  in  the  streamlines  plotted  in  Figure  133.  The  dividing  streamline  y/  =  0 
extends  approximately  0.5^ above  the  surface,  and  reattaches  at  x/S=\.9.  The  modification  of 
the  time-averaged  velocity  profiles  is  shown  in  Figure  IV.20  in  wall  coordinates.  (The  friction 
velocity  of  the  baseline  flow  at  the  actuator  orifice,  Ut=  0.405  m/s,  is  used  in  the  normalization). 
The  time-averaged  velocity  field  shares  many  features  with  the  quasi-steady  recirculation  region 
formed  at  a  much  higher  dimensionless  frequency  (§IV.3).  Within  the  interaction  domain 
{x/d=  1.0)  the  passage  of  the  discrete  vortices  induces  time-averaged  reversed  flow  adjacent  to 
the  wall.  The  presence  of  the  interaction  domain  displaces  the  cross-flow  fluid,  and  the 
associated  convective  acceleration  as  well  as  the  additional  momentum  injected  by  the  synthetic 
jet  results  in  a  maximum  in  the  profile  at  the  outer  edge  of  the  time  averaged  vortex. 
Downstream  from  the  mean  reattachment  point,  there  is  a  significant  wake  effect  in  the  velocity 
profile.  The  associated  velocity  deficit  reduces  as  a  new  inner  layer  within  the  boundary  layer  is 
formed  {x/d  =  4.0,  8.0).  In  fact,  the  velocity  distribution  at  x/S=  8.7  definitively  exhibits  a 
multi-tiered  character,  with  the  formation  of  an  inner  boundary  layer  (y^<100,  y/6<03) 
followed  by  a  region  of  constant  streamwise  velocity  (100  <  200,  0.3  <y/S<  0.6)  and  then 

an  outer  layer.  Under  these  conditions,  a  greater  distance  is  required  (>10^  for  the  velocity  field 
to  recover  a  normal  turbulent  boundary  layer  profile  compared  to  the  previous  case  (§IV.4). 

Within  the  interaction  domain,  all  three  components  of  Reynolds  stress  are  substantially  higher 
than  the  baseline  values  as  illustrated  in  Figure  IV.21.  Cross-stream  distributions  of  the  normal 
(mV  and  Vv')  and  shear  (mv')  components  of  Reynolds  stress  downstream  from  the  interaction 


62 


Figure  IVA9.  Left:  Phase-averagednormalizedvorticity  with  superimposed  velocity  vectors,  Re 0  =  700,  St S=  0.47,  Uj/U^ 

=  1.25,  S7b  -  6.5:  (a)  0°  (b)  60°  (c)  120°  (d)  180°  (e)  240°,and  ff)  300°  Right:  Time-averaged  normalized  vorticity  with 

velocity  vectors  (a  &b)  and  streamlines  (c  &d)  for  baseline  (a,  c)  and  actuatd  (b,d). 
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Figure  IK20.  Time-averaged  velocity  profiles  in  wall  units,  x/8  =  1. 0,2.0,  4. 0,  and  8. 7  (top  to  bottom):  baseline  (O),  actuated 
(%),  successive  locations  separated  by  one  decade. 


domain  are  shown  in  Figures  IV .22-Left  and  IV .22-Right,  respectively.  At  the  downstream  end 
of  the  interaction  domain  (x/S=^2.0,  Figure  IV.22-Left  a  and  IV.22-Right  a)  all  three  stress 
components  are  still  an  order  of  magnitude  above  the  baseline  values.  The  stresses  are  rapidly 
attenuated  with  downstream  distance  in  this  region  of  the  flow,  decreasing  by  a  factor  of  three  by 
x/S=  4.0  (Figure  IV.22-Left  b  and  IV.22-Right  b).  As  a  new  inner  region  develops  within  the 
boundary  layer  the  distributions  of  Reynolds  stresses  also  change,  as  local  peaks  develop  close  to 

the  wall  in  mV  and  17v'  atx/S=  8.7  (Figures  IV.22-Left  c  and  IV.22-Right  c).  The  present  data 
set  is  not  of  sufficient  extent  to  determine  the  streamwise  length  required  for  the  return  to  a 
normal  turbulent  boundary  layer. 

IV.5.2  The  Pressure  Field 

The  discrete  vortices  formed  by  the  interaction  between  a  synthetic  jet  and  cross  flow  induce 
large  pressure  gradients  locally.  The  normalized  phase-averaged  streamwise  <  dpfdx  >  and 
cross-stream  <  dp/dy  >  pressure  gradients  are  shown  for  six  instances  during  the  actuation  cycle 
for  the  current  experimental  conditions  in  Figure  IV.23.  As  observed  in  earlier  measurements, 
each  CW  vortex  produces  alternating  sign  gradients  (favorable,  then  adverse),  which  diminish 
rapidly  and  completely  disappear  within  two  wavelengths.  As  a  result,  there  is  a  finite  domain 
immediately  downstream  from  the  actuator  orifice  where  the  time-averaged  pressure  gradients 
are  modified  as  shown  in  Figure  IV.24  (computed  via  the  Reynolds-Averaged  Navier  Stokes 
equations).  The  time-averaged  pressure  gradients  bear  a  remarkable  resemblance  with  the  phase- 
averaged  gradients  associated  with  a  single  vortex,  suggesting  that  the  time-averaged  interaction 
may  be  modeled  as  a  stationary  vortex.  The  normalized  time-averaged  pressure  field,  shown  in 
Figure  IV.24,  is  determined  by  integrating  the  pressure  gradients  using  the  combined-path 
Runge-Kutta  routine  outlined  in  Appendix  B.  A  low-pressure  region  is  produced  immediately 
downstream  from  the  actuator  orifice,  with  a  minimum  pressure  of  p--2.\. 

IV.5.3  Boundary  Layer  Evolution 

As  discussed  in  §IV.5.1,  under  the  present  conditions  the  combined  action  of  the  passage  of  the 
discrete  vortices  results  in  the  formation  of  a  time-averaged  stationary  vortex,  which  is 
sufficiently  strong  to  produce  a  large  time-averaged  recirculation  region.  One  indication  of  the 
re-establishment  of  a  new  inner  region  of  the  boundary  layer  is  seen  in  the  wall  shear  stress, 
computed  from  least-square  error  polynomial  curve  fits  to  cross-stream  profiles  of  the 
streamwise  velocity  component.  The  normalized  wall  shear  stress  is  shown  at  six  phases  of 
actuation  in  Figure  IV.25  along  with  the  time-averaged  actuated  and  baseline  values.  In  the 
baseline  flow,  the  wall  shear  at  the  actuator  location  is  f„=  0.0033.  At  the  onset  of  blowing  (<^= 
0°),  there  is  a  local  minimum  in  at  x/S=  1.1  associated  with  the  previously  formed  CW 

vortex.  During  the  blowing  portion  of  the  cycle  (^=  60°,  120°)  the  wall  shear  stress  upstream 
from  the  actuator  orifice  becomes  smaller,  and  large  negative  shear  is  produced  as  the  vortex  pair 
rolls-up  (at  60°,  f„=  -0.012  atx/^=  0.16,  compared  with  the  baseline  value  of  0.0033).  As 
expected,  the  deviation  in  shear  stress  from  the  time-averaged  value  associated  with  each  vortex 
rapidly  diminishes  downstream.  During  the  suction  portion  of  the  cycle  (^  =  240°,  300°)  the 
shear  stress  is  upstream  of  the  orifice  increased,  and  decreased  immediately  downstream  as  fluid 
is  drawn  into  the  actuator.  The  net  effect  is  a  region  of  reduced  wall  shear  stress  immediately 
downstream  from  the  orifice,  in  fact,  negative  shear  persists  for  0  <  x/S^  1 .9  as  a  mean 
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Figure  IV.21.  Reynolds  Stress,  baseline:  uu'lul,  vVful  &  u'v'Jul  baseline  (a,  b  &c)  and  actuated  (d,  e,  &J). 


Figure  IV.22.  Cross-stream  distribution  ofJLsji:  uuful,  &  vv'lul :  (a)  x/S  =  2.0,  (b)  4.0,  and  (c)  8.67:  (baseline  in  solid  and 
dashed  lines,  and  actuated  •  &  M),  Right:  uV  Ul ,  x/S=  2.0  (a),  4.0  (b),  and  8.67  (c).  baseline  in  solid  line  and  actuated  in  A). 
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Figure  IV.23.  Phase-averaged  normalized  pressure  gradients,  C  (a-fi;  <  dpjdy  >  (g-l):  (a,  g)  ^  ~  0^,  (b,  h) 
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Figure  IV,24,  Normalized  Reynolds-averaged  pressure  gradients,  dpjdx  dpj^y  baseline  (a,  b)  and 
actuated:(c,  d).  Normalized  time-average  pressure,  baseline  (d)  and  actuated  (e)  (Dashed  and  solid  lines 
indicate  least  significant  contours,) 


66 


recirculation  region  is  created  (cf.  Figure  133).  Downstream  from  this  region,  the  shear  stress  in 
the  presence  of  actuation  rapidly  recovers  its  baseline  value  as  a  new  inner  region  in  the 
boundary  layer  is  formed. 

Both  the  displacement  thickness  (^)  and  momentum  thickness  {d)  are  substantially  modified 
across  the  entire  measurement  domain  by  the  interaction,  as  shown  in  Figure  IV.26.  The  discrete 
vortices  displace  the  cross  flow  fluid  within  the  interaction  domain,  causing  a  local  maximum  in 
S  at  x/S=  1.  Within  the  same  region,  the  additional  momentum  input  of  the  jet,  as  well  as  the 
convective  acceleration  of  the  cross-flow  due  to  the  displacement  effect  results  in  negative 
momentum  thickness  (indicating  that  the  flow  contains  more  momentum  than  a  uniform  profile 
at  C/oo).  The  shape  factor  hn,  which  is  has  little  physical  meaning  within  the  interaction  domain, 
is  reduced  downstream  fi’om  the  interaction  domain  compared  to  the  baseline  {hn  =  1.40  vs. 
1.55),  indicating  that  the  boundary  layer  has  more  momentum  immediately  adjacent  to  the  wall 
as  a  new  inner  layer  is  formed  (cf.  Figure  IV.20).  It  is  apparent  that  compared  to  the  case  with 
discrete  vortices  scaling  with  the  inner  portion  of  the  boundary  layer  (§IV.4)  the  boundary  layer 
characteristics  downstream  from  the  interaction  domain  require  a  much  greater  streamwise 
distance  to  fully  evolve. 

IV.5.4  Turbulent  Kinetic  Energy 

The  normalized  time  averaged  turbulent  kinetic  energy  (^^ )  is  plotted  alongside  the  normalized 

rates  of  turbulence  production  (P)  and  viscous  dissipation  (£)  in  Figure  IV.27  for  the  baseline 
and  actuated  flows.  In  the  baseline  flow,  turbulent  kinetic  energy  is  distributed  throughout  the 

boundary  layer  (profiles  of  uu  and  vV^are  shown  in  Figure  IV.2)  while  the  highest  levels  of 
production  and  viscous  dissipation  are  adjacent  to  the  wall.  In  the  actuated  flow,  the  largest 
magnitudes  of  ,  P ,  and  e  are  confined  to  the  interaction  domain  and  are  an  order  of 
magnitude  larger  than  those  present  in  the  baseline  flow.  Particularly  large  values  for  all  three 
quantities  occur  immediately  above  the  actuator  orifice,  where  the  velocity  fluctuations 
introduced  by  the  synthetic  jet  are  largest.  Negative  turbulence  production,  occurring  above  the 
actuator  orifice,  indicates  a  transfer  of  energy  from  the  periodic  jet  motion,  into  the  mean  flow. 
The  largest  [positive]  turbulence  production  occurs  in  the  domain  where  the  discrete  vortices 
break  down  (0.2  <  x/8  <  2).  Downstream  from  the  interaction  domain,  the  largest  rates  of 
production  and  dissipation  occur  adjacent  to  the  surface  within  the  newly  formed  inner  layer 
{y/5  <  0.2,  cf.  Figure  IV.20). 

Power  spectra  of  the  u  velocity  component  are  obtained  from  LDV  data.  The  streamwise 
variation  of  power  spectra  at  constant  elevation  within  the  boundary  layer  (y/8  =  0.32)  is  shown 
in  Figure  IV.28a,  for  the  baseline  and  actuated  flows.  Within  the  interaction  domain  (pc/S  =  0.8, 
1.5)  the  entire  spectra  is  shifted  since  there  is  substantially  more  turbulent  kinetic  energy  in  the 
flow.  In  addition  to  the  spectral  peak  at  the  actuation  frequency  (kS  =  2.9,  f  =  300  Hz),  higher 
harmonics  are  also  present  within  this  region.  At  x/S  =1.5  the  interaction  leads  to  an  expanded 
inertial  subrange,  with  a  -5/3  slope  starting  at  the  actuation  frequency  and  extending  at  least  a 
decade.  In  the  region  downstream  from  the  interaction  domain,  where  a  normal  turbulent 
boundary  layer  begins  to  redevelop,  the  spectra  contain  multiple  changes  in  slope  as  it 
redevelops  a  normal  energy  distribution.  While  the  spectral  peak  at  the  actuation  frequency 
decreases  with  downstream  distance  from  the  actuator,  it  is  still  present  at  x/S  =  IS. 
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Figure  IV.2S.  Normalized  wall  shear  stress,  ^  =  0°,  60“’.  120°,  180°,  240°,  and  300°  (from  top  to  bottom, 
successive  phases  separated  by  0.02):  (baseline-gray  line:  time-averaged  actuated-black  solid  line,  phase- 
averaged  A/ 


Figure  IV.26.  Boundary  IcQfer  evolution,  baseline:  S*  (O),  0(0),  hi 2  (A);  actuated:  S*  (%),  0(M),  hi 2  (A). 


Figure  IV.28.  u power  spectra:  a)y/6=  0.32  (y+  =  85),  x/S=  0.8,  1.5,  3.0,  6.1,  12,  and  18  (top  to  bottom),  b)x/S 
=  6.1,  y/8=  0.08,  0.16,  0.32,  and 0.64  (top  to  bottom)  fy+  =  21,  43,  85,  and  170),  c)x/5=  12,  y/S=  0.08,  0.16, 
0.32,  and 0.64  (top  to  bottom)  =21,  43,  85,  and  170),  andd)x/S=  18,y/S=0.08,  0.16,  0.32,  and  0.64  (top  to 
bottom)  (y+  =  21,  43,  85,  and  170.  Baseline  in  gray,  actuated  solid  line.  Successive  spectra  displaced  by  two 
decades. 


The  variation  in  u  power  spectra  across  the  boundary  layer  is  shown  at  x/S  —  6.\,  12,  and  18  in 
Figures  IV.28b,  lV.28c,  and  IV.28d,  respectively.  At  x/S  =  6.1,  the  actuated  flow  still  contains 
more  turbulent  kinetic  energy  than  the  baseline,  particularly  so  at  higher  cross-stream  elevations 
above  the  surface.  Curiously,  in  contrast  to  the  baseline  flow  the  power  spectral  density  of  the 
actuated  flow  is  roughly  constant  for  <  1,  and  has  approximately  a  -f  slope  at  higher 
frequencies.  Further  downstream  (x/S  =  12,  18),  except  for  a  small  peak  at  the  actuation 
frequency,  the  baseline  actuated  spectra  are  almost  indistinguishable  for  y/S  <  0.16  <  43), 

which  is  within  the  inner  tier  of  the  redeveloping  boundary  layer.  The  outer  part  of  the  boundary 
layer  requires  greater  streamwise  distance  to  recover  a  normal  turbulence  structure  due  to  the 
wake  effect  caused  by  the  presence  of  the  interaction  domain. 

IV.6  CONCLUSION:  DEMONSTRATION  OF  SEPARATION  CONTROL 

As  demonstrated  in  earlier  work  at  Georgia  Tech,  synthetic  jet  actuation  can  significantly  alter 
the  aerodynamic  characteristics  of  a  circular  cylinder.  In  particular,  the  modification  of  the 
streamwise  pressure  gradient  within  the  interaction  domain  leads  to  a  thinner  boundary  layer 
downstream,  which  is  more  capable  of  withstanding  adverse  pressure  gradients.  The  resultant 
delay  in  boundary  layer  separation  leads  to  the  generation  of  lift,  and  the  reduction  in  pressure 
drag.  In  this  section,  the  effects  of  synthetic  jet  actuation  on  a  separated  flow  are  studied  within 
the  flat  plate  boundary  layer. 

A  separated  flow  region  on  the  flat  plate  was  created  by  mounting  a  75  mm  tall  obstruction  with 
a  downstream  fairing  to  the  flat  plate  model.  The  free  stream  velocity  and  actuation  conditions 
are  identical  to  those  in  §IV.4  (C/«=  10  m/s,/=  1130  Hz,  Uj=  12.5  m/s,  and  6  =  0.051  mm). 

The  streamwise  location  (125  mm  downstream  from  the  orifice)  was  selected  such  that  boundary 
layer  separation  occurs  downstream  of  the  interaction  domain.  In  the  absence  of  actuation,  the 
wall-mounted  obstruction  induces  the  formation  of  a  large  separated  flow  region  and  a  trapped 
vortex,  as  seen  in  Figures  IV.29a  (time-averaged  velocity  distributions)  and  149a  (corresponding 
streamlines).  Jet  actuation  leads  to  the  formation  of  a  time-averaged  separation  bubble  that  is 
characterized  by  the  dividing  streamline  yf=  0  (Figure  149b).  The  presence  of  the  Interaction 
domain  displaces  and  accelerates  the  cross  flow  fluid  (Figure  IV.29b).  Remarkably,  even  in  the 
presence  of  such  a  large  adverse  pressure  gradient  (imposed  by  the  obstruction)  a  new  inner  layer 
is  formed  within  the  boundary  layer,  as  seen  in  the  velocity  distribution  at  x/h=  140.  The 
creation  of  this  inner  layer,  ostensibly  the  result  of  a  modification  of  the  streamwise  pressure 
gradient  within  the  interaction  domain,  results  in  a  downstream  movement  of  the  separation 
streamline  from  x/h  =  1 1 0  to  x/h  =  140. 
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